AD-A099  191  SYSTEMS  CONTROL  INC  <VT>  PALO  ALTO  CA  AERONAUTICAL  A— ETC 
PARAMETER  INSENSITIVE  CONTROL. (U) 

APR  81  R  L  KOSUT 


N00014-77-C 


uiiu  till  uurr  nu 


AERONAUTICAL  £>  MARINE  SaSTEMS 

DIVISION 


SsSTtMS  CONTROL,  INC.  (VO  ■  1801  PAGE  MILL  DOAD  ■  PALO  ALTO,  CA  9A30A  ■  TELEX  3A8-A33  ■  (AI5>  A9A-U65 


y  -x  ;y  ■  -/•  - '/ 


APRIL  1<981 

*  .... 


FINAL  REPORT 

~fC 


PARAMETER  INSENSITIVE  CONTROL  - 

■"  jjfr 

Project  5203 


Prepared  by 
Robert  Ly^Kosut j 


Prepared  for  the 
OFFICE  OF  NAVAL  RESEARCH 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


e>-#v 

G 


_ — -r  =;«•  >\  ■t  V  VJ’.TtT  ~ A  y 


,y  ter  ?■'  ■  ■  ' 


appioy'x1  > 
rn.’P  -ay. _ 


XtJ 


/  /  — 


Si  4  27  061 


TABLE  OF  CONTENTS 


Page 


I.  INTRODUCTION  AND  SUMMARY  .  i 

II.  ROBUSTNESS  OF  THE  LQG  CONTROL  .  7 

III.  MEASURES  UF  ROBUSTNESS  .  15 

3.1  Time-Domain  Measures  of  Kobustness  .  14 

3.1.1  Linearized  Control  .  lo 

3.1.2  Reduced -Order  Control  .  j.8 

3.2  Frequency- Doma in  Measures  of  Rooustness  .  la 

3.2.1  Unmodeled  Perturbations  .  2 

3.2.2  Reducea-Order  Control  .  2 

3.2.3  Actuator/Sensor  Dynamics  .  26 

REFERENCES  .  29 

APPENDIX  A:  LINEAR  PARAMETER  INSENSITIVE  CONTROLLER^  AND 

FILTERS  .  51 

APPENDIX  B:  FREQUENCY-SHAPED  COST  FUNCTIONALS:  AN 

EXTENSION  OF  LQG  DESIGN  METHODS  .  65 


APPENDIX  C:  THE  ACCOMMODATION  OF  LARGE  PLANT  PARAMETER 


VARIATIONS  BY  STATE  FEEDBACK  .  9" 

APPENDIX  D:  CONTROL  OF  FLEXIBLE  STRUCTURES  .  125 

D.l:  ROBUST  INVERSE  OPTIMAL  CONTROL  FOR  . 

FLEXIBLE  STRUCTURES  . ' .  125 

D . 2 :  STABILITY  OF  LQG  MODAL  CONTROL  FOR  . 

FLEXIBLE  STRUCTURES  .  14" 


hoc \ on  For 
nt’?  c"- 9 1 


i 

! 

I 

1 


I 


i  i  i 


CO  L/i 


I.  INTRODUCTION  AND  SUMMARY 


The  objective  of  this  research  study  is  to  understand  more 
about  the  characteristics  of  control  systems  that  are  parameter 
insensitive.  Consider  the  generic  feedback  control  system  as 
illustrated  in  Figure  1.  The  requirement  on  feedback  control  is 
to  maintain  system  performance  within  a  specified  tolerable  range 
irrespective  of  uncertainty  in  the  models  of  the  system  as  well 
as  the  disturbances.  The  reason  for  this  requirement  is  that  the 
control  is  desi gned  on  the  basis  of  a  nominal  model  of  the 
system,  whereas  the  control  operates  on  the  real  system  for  which 
the  model  is  only  an  approximation.  When  tne  true  system 
parameters  or  behavior  exceed  those  accounted  for  in  the  model 
the  controlled  system  may  fail  to  operate  within  the  required 
performance  envelope.  Thus,  the  need  for  parameter  insensitive 
control . 

The  generic  problem  described  above  gives  rise  to  two 
specific  problems  -  a  sensitivity  problem  and  a  robustness 
problem . 

The  sensitivity  problem  is  concerned  with  small  changes  in 
parameter  values  about  nominal  design  values,  and  the  relation  of 
these  small  changes  to  specific  system  properties  (.e.g. 
stability,  pole-placement,  etc.)  which  are  to  be  maintained 
within  specified  ranges. 

The  robustness  problem  is  concerned  with  classes  of  models 
or  sets  of  parameters  (including  the  nominal  class  or  set)  for 
which  specific  system  properties  (like  stability)  are  to  be 
maintained.  Clearly,  robustness  includes  sensitivity,  whereas 
sensitivity  does  not  imply  robustness. 

The  focus  in  this  study  is  on  robustness  of  linear  systems 
particularly  stability  robustness  and  performance  robustness. 
Stability  robustness  refers  to  the  ability  of  the  control  to 


maintain  staoilicy  in  the  presence  of  model  uncertainty. 
Performance  robustness  refers  to  tne  ability  of  tne  control  to 
maintain  a  measure  tor  measures j  of  performance  within  a 
specified  tolerable  range.  Stability  rooustness  nas  received 
considerable  attention,  for  example  see  the  discussion  in  Lij  anc 
the  detailed  coverage  in  l  3 j .  Performance  robustness  aas  not 
received  the  same  attention,  primarily  aue  to  the  general  lack  of 
satisfactory  results  in  stability  robustness,  ana  the 
considerably  more  difficult  nature  of  the  per formanc e  robustness 
problem.  Attempts  have  been  made  in  performance  robustness,  for 
example,  maintaining  a  prescribed  degree  of  exponential  staoility 
[3]  in  the  presence  of  uncertainty.  In  addition,  the  focus  has 
been  on  analysis  of  robustness  ii.e.  how  rooust  is  the  control) 
rather  than  on  the  syntnesis  of  robustness  li.e.  the  design  of 
control  with  a  specified  measure  of  robustness).  The  ultimate 
aim  would  be  a  method  to  synthesize  control  laws  wnicn  have 
specific  stability  robustness  and  performance  rooustness 
properties. 

New  approaches  to  the  stability  ana  performance  rooustness 
synthesis  problems  are  presented  in  this  report  as  well  as  areas 
for  future  research. 

The  appendices  contain  detailed  investigations  in  several 
areas.  Appendix  A  discusses  the  performance  robustness  issue  in 
a  very  precise  setting  by  defining  a  worst  case  design  ana  a  best 
average  design  optimization  procedure.  It  is  shown  that 
state-feedback  or  observer-based  feedback  controllers  with 
constant  gains  are  not  optimal  when  the  parameters  of  the  sytem 
change.  An  adaptive  control,  which  is  optimal,  updates  the  gains 
in  accordance  with  parameter  changes  by  effectively  using  past 
knowledge  of  the  states.  The  controllers  ana  filters  proposed 
here  use  a  limned  amount  of  past  states  but  with  constant 
gains.  The  gains  are  selected  to  either  maximize  worst  case 
performance  or  maximize  average  performance. 


Appendix  B  snows  new  results  in  cue  use  or  LQG  design  oy 
extending  the  theory  to  allow  tor  f requency-snaping  in  tne  oost 
functional.  This  extention  proviaes  a  means  to  meet  classical 
robustness  requirements  with  the  automated  computational 
procedure  of  modern  control  theory.  Tne  f requency-snaping  metnou 
permits  the  designer  to  synthesize  an  optimal  compensator  winch 
takes  into  account  a  priori  Knowledge  of  disturbances  and 
parameter  uncertainty.  For  example,  if  the  transfer  function  or 
the  model  weie  known  to  lie  within  a  certain  range  of  transfer 
functions  (corresponding  to  a  range  of  parameter  variations!  tnen 
the  cost  functional  can  be  frequency  snapeo  so  that  equal  penalty 
is  applied  across  this  range.  The  degree  of  penalty  must  oe 
determined  iteratively  by  evaluating  the  performance.  This 
procedure  produces  performance  robustness  to  a  specified  range  or 
parameter  variations  or  disturoances .  Appendix  B  snows  four 
examples  of  applying  this  methodology:  flexible  spacecraft, 
helicopter  vibration,  industrial  crane,  and  aircrart  in  lateral 
wind . 

Appendix  C  approaches  the  robustness  proolem  from  an 
algebraic-geometric  view.  Stable  suospaces  are  designed  for  a 
class  of  controls  whicn  remain  stable  for  a  specified  class  of 
parameter  variations.  This  analysis  also  discriminates  systems 
where  linear  state  feedback  control  laws  work  adequately  and 
systems  where  adaptive  control  laws  are  required. 

Appendix  D  presents  methods  to  design  robust  controllers  for 
flexible  structures.  Flexible  structures  such  as  an  aircraft 
wing,  long  beam,  or  a  large  space  structure  (LSd)  are 
representatives  from  a  class  of  systems  with  infinite  mooes.  Tne 
actual  dynamics  are  governed  by  partial  differential  equations 
for  which  the  lower  frequency  modes  can  be  adequately  (in  most 
cases)  modeled  by  interconnected  oscillatory  systems.  The 
difficulty  arises  in  that  the  dimension  ot  this  model  is 
extremely  large  (e.g.,  thousands  of  inodes!,  and  the  higner  modal 
frequencies  and  modes  shapes  are  poorly  known.  However,  tne 


4 


structure  ot  the  >nouel  is  definitely  known--!. e.,  positive 
symmetric  matrix  second  oruer  equations.  The  prooiem  then  is  to 
find  robust  controllers  whicn  are  based  on  a  lower  order  mouei  ot 
the  flexible  structure.  Two  approaches  are  presented.  The  first 
method  (Appendix  D.l)  uses  the  solution  of  an  inverse  optimal 
controller  problem  to  estaolish  an  iterative  procedure  to  design 
robust  co-locatea  vibration  control,  for  pooriy  uampea  tlexiole 
structures.  The  solution  requires  output  feeaoack  witn  specnieu 
constraints,  leading  to  rooustness  with  respect  to  unmoueled 
moues  and  a  large  class  of  parameter  variations.  The  rooustness 
properties  are  proved  directly  from  Known  properties  of  control 
laws  optimising  quadratic  perrormance  measures.  Tne  second 
method  (Appendix  D-2J  uses  the  frequency  domain  rooustness 
measures  discussed  in  Section  3.2,  to  design  an  LqG  modal  control 
which  is  robust  with  respect  to  tne  unmoueled  but  bounded 
residual  modes  of  the  structure.  The  rooustness  measures  are 
used  to  synthesize  a  snape/vioration  LQG  moaal  control  using 
position  measurements. 

1.1  SHORT  SURVEY  AND  FURTHER  ISSUES 

Section  II  describes  some  of  the  issues  of  robustness  with 
respect  to  the  well  established  LQG  control.  Section  III  surveys 
recently  developed  measures  of  robustness  in  the  time-domain  and 
in  the  frequency-domain.  Several  promising  areas  for  future 
application  and  investigation  of  these  robustness  metnous  are 
also  presented  which  hignlignt  new  ways  to  consider  the  issues  in 
linearized  gain  scheduling  control,  reduced  order  control,  the 
effect  of  unmodeled  but  bounded  perturbations,  and  tne  effect  of 
uncertainty  and  reliability  of  actuators  ano  sensors. 
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II.  ROBUSTNESS  OF  THE  LQG  CONTROL 


Consider  the  system  governed  by  the  t ime- invar i ant  linear 
model , 


x  =  Ax  +  Bu  +  w  ,  t^_o 
y  =  Cx  +  v 


(1) 


where  x  is  the  state,  u  the  control,  y  the  sensed  output,  and 
(w,  v)  the  dynamic  and  measurement  disturbances,  respectively. 

The  disturbances  (w,  v)  are  usually  assumed  to  be  zero-mean, 
white,  Gaussian  processes  with  known  constant  intensities.  For 
purposes  here  it  is  sufficient  to  assume  simply  that  they  are 
bounded*  and  may  have  known  spectra.  Since  y  is  the  available 
measurement,  the  control  is  restricted  to  the  form, 


u  =  -Hcy 


(2) 


where  Hc  is  a  linear,  t ime- invar iant  operator  with  rational 
transfer  function  Hc(s).  The  feedback  structure  is  shown  in 
Figure  2. 

It  is  well  established  [4]  that  the  LQG  control, 
u  =  -  Fx 

x  =  Ax  +  Bu  +  K(y-Cx)  (3) 

stabilizes  the  system  provided  that  (v,  w)  are  bounded  and  the 
(F,  K)  gain  matrices  are  selected  so  that  the  eigenvalues  of 
both  A-BF  and  A-KC  all  have  negative  real  parts.  If  (A,  B) 
is  controllable  and  (A,  C)  is  observable  then  the  eigenvalues 
of  A-BF  and  A-KC  can  be  arbitrarily  assigned  [4J.  Stability  is 
also  guaranteed  if  (F,  K)  are  computed  from 

F  =  BTP  ,  K  =  SCT 

*  the  expression  Hx  Ts  bounded"  as  used  here  means, 


/„ 


xl(t)  x  (t)  dt  <  =,  i.  e.  x  is  in 


L,  unless 


otherwise  stated. 
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where  P  and  5  are  the  symetr ic  positive  definite  solutions  it 
tne  algebraic  R i c a 1 1 i  equations. 
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The  source  of  the  robustness  problems  of  LQG  control  can  be 
simply  stated:  There  is  no  guarantee  tnat  the  eigenvalues  of  the 

LOG  comoensator  are  stable. 

- - - * -  - 

Reference  5  gives  an  excellent  description  of  this  prooiem. 
The  results  of  instability  of  the  compensator  H  can  be 
disastrous.  Consider  what  happens  if  a  control  loop  or  a  sensor 
becomes  disabled:  The  system  will  most  likely  be  unstable. 

Also,  changes  in  the  true  system  model  (A,  B,  C)  matrices  may- 
cause  instability.  This  results  because  no  measure  of  robustness 
has  been  intentionally  designed  into  the  system. 

However,  it  has  been  shown  l 4 , 6 , 7 J  that  tne  full  state 
feedback  LQ  control, 

u  =  -Fx  l 9) 

with  F  from  (4)-(5),  guarantees  stability  for  1/2  to  infinite 
gain  margin  per  channel,  or  +  60°  pnase  margin  per  channel. 
Unfortunately,  this  guarantee  is  too  restrictive  in  regard  to 
robustness  and  other  practical  considerations  of  multi-loop 
control.  First  of  all,  LQ  control  requires  full  state  feedback 
into  all  input  channels.  Generally  this  is  impossible,  since  ail 
the  states  are  not  available  for  measurement.  Secondly,  the 
closed-loop  system  with  LQ  control  rolls  off  at  the  high 
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trequencies  at  6  do/octave.  Thus,  any  unmoueieu  high  frequency 
phenomena  in  the  system  can  have  a  deleterious  ettect  on 
performance  and  stability.  Tmruiy,  ana  most  important  to 
noustness,  most  parameter  variations  uo  not  appear  in  tne  torn 
of  gain  or  pnase  perturoat  ions .  Thus,  the  gain  arid  phase  margin 
properties,  attributable  to  tne  Lq  control,  are  not  vanu 
measures  of  rooustness. 

.An  even  more  distressing  situation  arises  wnen  utilizing 
observer -based  LQG  control  (5J.  In  this  case,  there  are  no 
guaranteed  margins  at  alii  However,  the  observer  gain  K  can  oe 
adjusted  [81,  under  certain  conditions,  so  tnat  the  gam  ana 
pnase  margins  of  the  LQG  control  approach  asymptotically  the  gain 
and  phase  margin  properties  of  the  LQ  control.  But  this  does  not 
mitigate  the  robustness  proolem  since  gain/phase  margins,  as 
already  mentioned,  are  not  necessarily  measures  of  rooustness. 

However,  the  gain/phase  margin  properties  of  the  LQ  control 
result  from  a  more  general  consideration  discussed  in  [7J. 
Consider  the  perturbed  LQ  control  shown  in  Figure  5.  The 
perturbation  L  is  assumed  to  oe  a  bounaea  linear, 
t ime - invar iant ,  operator  with  rational  transfer  function  Lisj. 
Under  these  conditions  if 

L(j)+LT(-j)  I,  0  [iUJ 

then  the  perturbed  LQ  control  is  stable.  The  gain  ana  phase 
margin  properties  of  LQ  control  can  oe  derived  from  [10). 

However,  (10}  allows  for  a  broader  class  of  uncertainties,  for 
example,  perturbations  in  the  B  matrix  of  (lj.  For  perturbations 
in  A  and  C  of  (1J,  a  more  general  measure  of  robustness  is 
r e  qu i red . 
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III.  MEASURES  OF  ROBUSTNESS 


Perhaps  the  best  and  most  inclusive  measures  of  robustness 
are  those  pioneered  by  Nyquist  [9]  and  Bode  [10J  for  the 
s ing le - input ,  s ing le -output  control  system.  For  example,  the 
effects  of  plant  parameter  variations  or  unmodeled  high  frequency 
phenomena  are  immediately  apparent  in  the  frequency  domain  using 
Nyquist  plots.  The  gain  and  phase  margins  can  be  quantitatively 
established  by  inspection.  If  these  margins  are  not  adequate 
with  a  given  control  structure,  then  various  combinations  of 
lead/lag  compensation  networks  can  be  usea  to  adjust  the 
margins.  However,  as  soon  as  there  is  more  than  one  input/output 
pair,  the  simple  graphic  understanding  of  the  Nyquist  plot 
vanishes  in  a  puff  of  dimensionality,  and  the  problem  of 
robustness  becomes  more  difficult.  As  was  seen  by  the  previous 
discussion  of  LQG  control,  the  robustness  problem  remains,  even 
with  the  move  away  from  the  frequency  domain  and  the  introduction 
of  state-space  and  optimal  control  by  Kalman  [11].  Of  course  the 
1  inear -opt ima 1  time-domain  approach  does  yield  a  stable  control 
for  most  multivariable  systems,  provided  relatively  mild 
conditions  are  satisfied  (like  controllability  and 
observab i 1 i ty .  ) 

Multivariable  control  can  also  be  designed  in  the 
frequency-domain.  The  approaches  in  references  [12],  [15] 
require  the  concept  of  diagonal  dominance.  The  control  is 
designed  in  two  parts.  First,  the  system  is  diagnalized 
(abstractly)  and  secondly  the  "classical"  control  concepts  are 
applied  to  the  independent  diagonal  loops.  This  is  a  cumbersome 
approach,  but  even  so,  there  is  no  guarantee  of  robustness  since 
the  actual  control  is  not  diagonal  (there  are  transformations 
into  and  out  of  the  "diagonal"  space). 
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The  choice  of  designing  the  control  in  the  frequency-domain 
(the  so  called  "classical"  approach)  or  in  the  tiine-aomain  (the 
so  called  "modern"  approach)  is  all  wrapped  up  with  the  question 
of  robustness.  The  effect  of  uncertainty,  which  the  designer 
knows  is  "out  there"  may  be  more  describable  in  either  the 
time-domain  or  the  frequency-domain.  Attempts  to  measure  and 
test  robustness  have  been  proposed  for  both  the  time-domain  and 
the  frequency-domain. 

3.1  TIME- DOMAIN  MEASURES  OF  ROBUSTNESS 

Methods  have  been  developed  to  measure  staoility  robustness 
in  the  time-domain  (14].  Consider  the  system, 

x  =  (H  +  o H)x  (10) 

where  H  is  stable  and  6H  is  a  perturbation  in  the  model.  The 
system  (10)  is  stability  robust  if  and  only  if  the  eigenvalues 
of  H  +  3H  all  have  negative  real  parts.  This  requires  specific 
knowledge  of  5 H  which  may  not  be  available.  For  example, 
only  j  j 6H  ||  may  be  known.  The  approach  in  114]  is  to  form  the 


Lyapunov 

function. 

< 

II 

X 

H 

T) 

X 

(11) 

where  P 

-  P^>0  and  Q  =  >0  satisfy, 

and 

HTP  +  PH  +  Q  =  0 

V  =  -xTQ  x  +  2  xTP  Hx 

i-l  i  x|  l2(  -  X  (Q)  +  2X(p)  |  |  5  H  I  !  ) 

(12) 

is  always  negative  if, 

!  I  <5  H  j  |  <  (13) 

2~xJD 


Consequently  if  (13)  is  satisfied  then  V  is  a  Lyapunov  function 
and  so  (10)  is  asymptotically  stable.  (The  notation x  (• )  and 
X(»)  represent  the  minimum  and  maximum  eigenvalues, 
respectively,  of  the  matrix  argument.) 
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Consider  the  perturbed  system  of  (ij, 


x  =  (A  +  5  A )  x  +  (B  +  -5  B )  u  +  w 


y=(C+3C)x+v 
with  LQG  control  (3), 
u  =  -  F  x 

x  =  Ax  +  Bu  +  K(y-Cx) 

The  matrices  in  (10)  are  then, 

H  =  |a-bf  BF 
0  A-KC 

j  u  3A-5BF  5BF 

5A-K5C-5BF  5BF 

The  robustness  test  (13)  can  now  be  applied  using  the  above 
(15)  quantities.  Clearly  (15)  can  be  used  to  test  a  wide  variety 
of  parameter  variations,  like  those  in  (14).  There  are 
limitations.  For  example: 


U4  j 


(15) 


(1)  The  state  representation  restricts  the  form  of  the 
per turbat ion  (e . g . ,  unmodeled  high  frequency 
phenomena  can  not  be  accounted  for  in  (14j, 
without  extending  the  state  order.) 

(2)  If  the  state  order  is  too  large,  the  Lyapunov 
equation  (12)  may  not  be  solvable  numerically. 

(3)  It  is  not  clear  how  to  modify  the  control  if  the 
test  fails.  In  fact,  failure  of  the  test  does  not 
imply  instability.  The  test  (13)  is  sufficient 
for  stability. 


A  similar  test  exists  for  systems  with  non-linear  perturbations, 
x  =  Hx  +  h(x)  C 1 6 ) 


where  now, 

ii  h  ( x )  1 1  <  2l[.gj  llxll  (17) 

insures  asymptotic  stability  of  (16). 

Tests  (13)  or  (17)  can  also  provide  a  measure  of  performance 
robustness.  If  (P,Q)  satisfy, 

(H  +  a  I )  TP  +  P(H  +  oiI)+Q  =  0 

15 


with  a  >o,  then  (15)  or  (17)  implies  that  (lo)  or  vlbj  is 
guaranteed  to  be  stable  with  aegree  ot  exponential  stability 
equal  to  *  .  In  other  words,  the  nominal  system  witn  no 
perturbation  exhibits  performance  as  measured  by  the  j  -degree  ot 
exponential  stability,  and  the  system  will  be  performance  robust 
(in  the  sense  of  exponential  stability)  if  trie  perturoat  ions 
satisfy  (13)  or  (17)  . 

3.1.1  Linearized  Control 


An  interesting  potential  application  of  (17)  arises  from  the 
use  of  "linearized"  control  law  mechanizations.  A  common  design 
procedure  (e.g.  flight  control)  to  control  non-linear  systems  is 
to  linearize  about  several  equilibrium  conditions,  develop  tne 
linear  control  for  each  operating  point,  and  then  during  actual 
operation  schedule  the  appropriate  linear  gains.  The  real  issue 
with  this  approach  is  in  determining  a  minimum  number  of 
equilibrium  conditions  in  the  dynamic  envelope  and  a  gain 
scheduling  procedure  which  maintains  a  specified  degree  of 
performance  and  stability  throughout  the  dynamic  operational 
envelope.  Consider  the  non-linear  dynamic  system, 

x  =  f ( x , u) 

Define  Cx,  u  )  as  an  equilibrium  point  such  that 
o=  f(X,  u) 

In  general  there  are  an  infinite  number  of  such  equilibria. 
Suppose  one  candidate  equilibrium  has  been  selected.  Define  the 
errors , 

<5x  =  x  -  x 
6u  =  u  -  u 

consequently,  by  linearization  about  (x,  u), 
x  =  Ax  +  Bu  +  e 
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wue  re 


x  =  x 
u  =ii 


Select  the  control 
•u  =  -F 


B  = 


.  u 


X  =  X 

u  =  iT 


where  A-BF  is  stable.  A  fundamental  praolem  is  to  determine  a 
bound  on  '  x  to  insure  staoiiitv  or 

\  =  1  A - B F ) ■ x  *  e 


Applying  i  1  ■*  J  gives, 


where 


ana 


e 


_1Q) 
2  TIP) 


X 


(A-BF ) FP  *  P( A-BF )  ♦  0  =  0 
e  =  f  ( 7  ♦  :x,  u"  -  Fix)  -  (A-BF)  x 


The  problem  remains  to  evaluate  e  .  For  many  non-linear 
systems,  the  non-linearities  are  at  most  2nd  degree.  Thus, 

e  <  X  (  F  ,  X  ,  U  )  ;  :£X 

where  *  (  •)  is  the  norm  of  the  2nd  order  tensor  in  the  Taylor 
series  expansion  of  F(x  +  tx,  u  -  Fxx).  Consequently 

:i5x::<  rTTF7"ic!qo)-rp) 

determines  the  range  of  stability  of  the  ”1  inear i tea"  control. 
The  above  bound  is  only  illustrative  of  now  to  use  the 
time-domain  robustness  measure  (17).  The  complete  problem 
involves  many  equilibria,  gain  scheduling,  and  an  observer  since 
all  the  states  will  not  be  available  for  measurement. 
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5.1.2  Reduced  Order  Control 


.Another  potential  application  of  (131  is  to  ietermine  the 
stability  of  a  reduced  order  control.  Consider  the  controllable 
system, 

x  =  Ax  ♦  Bu 
u  =  -  Fx 

where  A-BF  is  stable.  It  dim(x)  is  very  large  then  a  reduced 
order  control  is  desired, 

u  =  -  k  ; 

:  =  Mx 

where  dimtc)  uim(x).  Usually  M  is  selected  by  insight  into 

the  svstem  behavior.  Generally  x  can  be  partitioned  so  that  x 
=  q).  In  this  case  K  is  just  those  gains  in  F  that 

multiply  the  truncated  states  : .  The  closed  loop  system  is, 

x  =  l  A  -  BK.'-i )  x 

Adding  ana  subtracting  BFx  gives, 
i  =  (A-BF)x  ♦  B( F-kMJx 
which  following  (13)  is  stable  it, 

BlF-kM;  <  TTT^T 

where  (P,  Q)  satisfy  (12)  with  H  =  A-BF. 

3.2  FREQUENCY-DOMAIN  MEASURES  OF  ROBUSTNESS 


Consider  the  stable  linear  unity  feedback  system  of  Figure  4 
Solving  for  e  gives, 

e  =  (I  +  GQ  )  * 1  u 

If  the  system  is  perturbed  as  in  Figure  5, 
then 


e  =  (  I 


3G) _1u 


J  and 


From  115],  [16],  and  [17]  it  u  is  bounded  (i.e.  u 

1 1  +  Gq  +  5  G) *  ^  is  bounded  then  the  perturbed  system 
is  stable.  But, 

(I  +  Gq  +  5G)'1  =  (I  +  II  +  Gq ) ' 1 i G } " 1  (I  ♦  G^'1 

Since  (I  +  G  )"^  is  bounded  (nominal  system  is  stable)  then 
stability  is  insured  if 

!  i  (I  +  (I  +  GQ  } ' 1 5  G )  * 1  '  !<  - 

Using  the  contraction  principal  [18]  the  aDove  is  bounded  if 

!  !  (  I  +  Gq ) ‘ 1 5  G  i  !  <  1  (18) 

This  is  proved  rigorously  in  [19].  Note  that  118)  uoes  not 
require  GQ  or  5G  to  be  linear,  t ime- invariant  operators. 

In  general, the  perturbations  can  be  classified  into  additive 
and  multiplicative  perturbations  as  shown  in  Figure  6. 

For  additive  perturbations  SG  =  L  and 

!  i  ( I  +  Gq ) *  1 L  !  !  <  1  (19a) 

insures  stability.  Isolating  the  additive  perturbation  gives  the 
more  conservative  condition, 

! !  L  j  |<  1/||(I+  Gq ) " 1  I  |  (19b) 

For  multiplicative  perturbations  5G  =  GQ  L  and 

!|(I  ♦  G  “1)“1  l  il<  1  (2  0a) 

o 

insures  stability.  Isolating  the  multiplicative  perturbation 
gives  the  more  conservative  condition, 

Ml  !]  <  1/ ! !  (i  ♦  go  " 1 )  ‘ 1  i  i  (zob) 
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It  Gq  and  L  are  linear,  t ime- i nva r ian t  operators  with  the 
cor  respond  ing  rational  transfer  functions  C>0  i  s  J  and  Lis) 
then  the  perturbed  system  rema ins  stable  if,  for  auditive 
per  tur bat  ions , 


'(  LCjuO  )<  c(  I  +  GQ  (.  jw  )  ), 


>  0 


121) 


or  for  multiplicative  perturbations, 

7  (  L(jw))c  £  (  I  +  G0_1(  joj))  , 


122) 


provided  that  the  perturbations  are  staole.  Relations  (21)  and 
(22)  are  proved  rigorously  in  [19]  and  also  in  (20)  where  these 
ideas  were  originally  developed  using  Nyquist  criteria 
arguments.  Note  that  o'  ( * )  and  o(*)  represent  the  maximum  ana 
minimum  singular  values  of  the  matrix  argument. 

Advantages  of  these  frequency-domain  robustness  measures 
over  their  time-domain  counter  parts  are  that  the  matrix  order 
of  Gq  is  much  lower  then  the  state  order,  and  unmodeled  high 
frequency  phenomena  are  easily  incorporated.  Tne  disaavantage  is 
that  G0  and  L  are  very  complicated  functions  of  a  complex 
variable . 


These  tests  can  be  applied  to  the  LQ  control.  For  the 
s ingle- input ,  single-output  case  Kalman  [ 6  J  showed  that  the 
magnitude  of  the  return-difference  for  an  optimal  control  is 

!  1  +  g0  C  )  !  £  1  ,  to  >  o 

Similarly  for  the  multi-loop  system  Anderson  and  Moore  l 4 ]  give 
the  more  general  result, 

(  I  +  Gq  ( -  ))  T(  I  +  G0  ( j  ^  I 

Both  imply  that 

a  (I  +  Gq)  >  1 

which  further  implies  that 
ad  +  G0'1)  >  1/2 

Consequently,  LQ  control  can  withstand  additive  perturbations, 

7  ( L)  <  1  ,  w  >  0 
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or  multiple  perturbations, 

T  Cl)  <1/2 

and  still  remain  stable.  In  general  these  results  uo  not  apply 
for  the  LOG  control. 

3.2.1  Unmodeled  Perturbations 

Consider  the  system  shown  in  Figure  7  where 

y  =  tHP  +  hr)  « 

u  =  Hc(v-y) 

!  v  i  i  <  =“ 

Hp  is  the  model  of  the  plant  dynamics  (.sometimes  referred  to 
as  the  primary  or  nominal  system  mouel); 

HR  is  the  residual  model,  where  usually  only  HRijui)|'  is 
known  (e.g.  high  frequency);  and  Hc  is  the  control  such  that 
the  loop  is  stable  with  HR  =  0.  Breaking  the  loop  at  u  gives, 

G  =  Hc  Hp  +  Hc  Hr 

Identifying  HqHr  as  an  additive  perturbation  gives 

a  (HCHR)<  2.  C I  +  HcHp),  u  >  0  125) 

as  a  sufficient  condition  for  stability.  If  Hc  represented 
an  LQG  design  then  Hc  would  have  to  be  stable  for  the  test  to 
work.  If  all  the  matrices  were  square  then  an  alternate  test 
results,  since 

G  =  HcHp(I  +  Hp'1  Hr) 

In  this  case  Hp  1  HR  can  be  identified  as  a  multiplicative 
perturbation  and  so, 

cCH'1  HB)  <  £(I  +  (HH)'1)  , -  >0  C  24 ) 

P  K  C  p  — 

insures  stability. 


in  1  -  0 


is  st a  Die  then 


t  il  e 


^  6  1  d  U  J  i  S  V  S  t  Cla 


/  it  n ,  i  s  stable  then  t  n  e  r  e  s  1  a  u  a  i  s  y  s  t  e;a 
tan  be  isolated  to  give  the  :n  ere  conservative  tone  it  ion, 

1 1 


l  ii^  ; 


H  ii  J 

C  D 


■) 


i  H .  J 


In  (.24)  tne  residual  system  ii,,  can  oe  isolated  ;  independent  ;f 
H„  stable.)  to  give  the  more  conservative  condition, 

■  iHn  )  -  l.H  )  •  ( I  + 1 H .  H  )  ^  .  •  0  ■  2  b  j 

•”  —  p  _  P  - 

Actual  computations  will  Determine  whicn  rooustness  measure 
of  the  above  are  the  least  conservative.  Of  course  i.  2  t> ;  can  oni 

be  used  if  H  ~  ^  exists  and  ( 2  5 can  onlv  be  usee  if  H 

p  '  w 

is  stable. 

Many  systems  fit  the  general  framework  o:  Figure 
Flexible  systems  (like  large  space  structures;  have  exactly  this 
format.  They  are  infinite  dimensional  systems  modeled  by  finite 
dimensional  svstems  with  much  iower  order  controllers.  In  this 
case  Hr  represents  the  effect  of  the  unmoaeieu  or  residual 
modes,  which  are  not  necessarily  ail  high,  frequency  effects.  T'n 
robustness  stability  of  LQG  modal  control  for  systems  of  this 
type  is  examined  in  detail  in  Appendix  D.2  using  the  frequency 
domain  robustness  measures  discussed  here. 


3.2.2  Reduced  Order  Control 


The  stability  of  the  reduced  order  control  discussed  m 
Section  1.3.2  can  also  be  examined  using  frequency  uomain 
robustness  measures.  For  that  system, 

-1 


G  =  KM  d  B  , 


(s I-A) 


Adding  and  subtracting 


v  -  p 


gives  , 


G  =  F  :  o  -  v  F  -  \M  )  :  3 

and  so 

ell  r -  KM  >  .'  3  i  l ,  i  *  r  :  B  J  ,  _  •  u 
insures  stabiiitv  or  the  reduced  order  contra 


i.d.3  Ac tua tor /Sensor  Dynamics 


Frequency -domain  robustness  measures  can  oe  used  to  examine 


the 

effect 

or  uncertainty  in  actuator/ sensor 

a  y  r* 

am  i  c  s  . 

dons  id 

the 

system 

of  Figure  8,  where  Hp  is  the  plant. 

H . 

A 

tne 

ac  t 

u  a  t  o  r  d 

ynamics,  the  sensor  dynamics, 

3  Pi  u 

u 

“C 

tne 

control. 

Suppose  that  is  perturoed  by 

'"‘A 

and 

repr  e 

sent  the  nominal  ac tua to r / s ensor  d> 

’  H  a. Ft 

1  cs  . 

dressing 

the  loop  at  the  actuator  input  gives, 

G  =  HcHsHpHA(I  *  Ha'X  i  Ha) 

Identifying  H^1  ;  HA  as  a  multiplicative  perturoat  ion  gives  the 
perturbation  bounds, 

T  (  '•  Ha)  v  lriA)  ^  II  +  ( 1  ;  ,  -  ^0 

■» 

which  insures  robust  stability,  provided  that  h.'L  :H,  is  staoie, 

ana  Hc  stabilizes  the  nominal  system.  A  similar  result  can 
also  be  established  for  sensor  perturbations. 

This  result  has  some  interesting  consequences.  Suppose  that 
for  the  nominal  system,  no  other  actuator/sensor  or  compensator 
gives  better  performance.  In  this  case,  the  robustness  test  can 
evaluate  candidate  actuators  which  are  not  optimal.  Tradeoffs  in 
actuator  parameters  and  robustness  are  easily  estaolished  without 
the  need  for  recalculating  a  new  optimal  control. 
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The  effect  of  actuator /sensor  parameter  uncertainty  and 
actuator/sensor  failure  can  also  be  tested.  If,  uncertainty 
results  from  actuator/sensor  failure  the  tests  give  an  indication 
of  reliability  to  failures  in  control/communication  channels  that 
are  not  catastrophic,  i.e.,  do  not  destabilize. 

It  should  be  pointed  out  that  where  the  loop  is  oroken  to 
calculate  G  is  dependent  on  the  uncertainty  being  tested.  A 
common  sense  rule  is  to  break  the  loop  wherever  nature  is 
inserting  a  disturbance  or  uncertainty.  There  is  no  sense 
testing  at  a  point  where  there  is  no  effective  disturbance. 
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APPENDIX  A:  LINEAR  PARAMETER  INSENSITIVE  CONTROLLERS  AND  FILTERS 


Several  methods  are  proposed  for  the  control  of  systems 
whose  mcdel  structure  or  model  parameters  are  not  known  exactly. 
One  class  of  control  laws  includes  explicit  real-time  estimation 
of  the  system  model  and  the  control  inputs  may  work  explicitly 
towards  enhancing  future  estimation  accuracy  (dual  control). 

These  control  designs  are  usually  complex  requiring  nonlinear 
controls,  even  for  linear  systems.  This  section  presents  a  broad 
class  of  linear  control  laws  which  are  minimally  sensitive  to 
variations  in  system  parameters.  It  is  shown  that  state-feedback 
controllers  with  constant  gains  are  not  optimal  when  the 
parameters  of  the  system  change.  Parameter  insensitive 
controllers  use  past  states  and  controls  for  feedback  and  on 
periodic  variations  in  control  gains.  These  controllers  do  not 
estimate  parameters  or  model  structure  and  the  system  need  not  be 
excited  to  improve  identif iability .  They  are  usually  simpler  to 
design  and  implement  than  adaptive  or  dual  controllers  and  often 
work  better  than  these  controllers  if  the  parameters  change 
relatively  quickly.  Classical  control  laws  are  shown  to  be  a 
subset  of  these  parameter  insensitive  controllers.  The  concepts 
and  formulations  of  parameter  insensitive  controllers  are  also 
extended  to  parameter  insensitive  filters.  Parameter  insensitive 
compensators  may  be  designed  by  combining  a  parameter  insensitive 
filter  and  a  parameter  insensitive  controller. 
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1.0  INTRODUCTION 


Consider  a  linear  system  in  which  Che  nxl  state  variable  x  follows 
the  difference  equacion 

x(k-t-l)  *  Ax(k)  •+■  But,k)  f  v(k;  '[lj 

u(k)  is  che  qxl  veccor  of  concrols  and  w(k)  is  an  nxl  vector  of  zero 
mean  which  Gaussian  state  disturbance.  A  and  3  are  state  transition 
and  control  distribution  matrices  of  appropriate  dimensions.  The  design 
of  control  inputs  for  this  system  and  its  continuous  time  counter  has 
been  a  problem  of  intense  study.  The  control  design  problem  consists  of 
selecting  the  control  input  sequence  u(k)  to  satisfy  certain  constraints 
on  states  and  controls  and  to  optimize  and  objective  function. 

Most  feedback  control  design  studies  have  been  limited  to  the  assumption 
that  A  and  3  are  known  perfectly  in  the  design  stage.  If  the  para¬ 
meters  of  the  system  are  known  exactly  and  ail  states  are  measured  with  no 
noise,  the  control  design  problem  simplifies  considerably.  In  addition, 
for  a  criterion  function  which  is  a  quadratic  function  of  state  and  control 
variables 


J  =  I  yT(k)Qy(k)  (2) 

k=*l 

y(k)  =■  Hx(k)  +  Du(k)  .  (3) 

the  optimal  control  input  at  any  time  is  a  linear  function  of  state  variables 
at  that  time.  This  control  law  has  several  desirable  features  including  stability 
or  the  closed-loop  system  if  (A,B)  and  (A . Q  )  are  controllable. 
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In  many  problems  of  practical  inCeresc,  all  model  parameters,  and 
sometimes  even  the  model  order,  are  not  known  exactly.  The  state  equation 
description  matrices  A  and  S  may  be  considered  functions  of  a  vector 
of  unknown  parameters  v,  whose  values  are  not  known  exactly  ac  the  time 
of  design  but  is  known  to  lie  in  a  set  .1.  3v  selecting  appropriately 
this  formulation  may  include  cases  where  the  model  order  is  not  known.  The 
control  design  problem  is  much  more  complex  in  this  case,  even  for  quadratic 
performance  indices.  Past  measurements  of  the  state  variable  are  then 
useful  for  controlling  the  system  in  the  future,  because  pasc  states 
contain  information  about  system  parameters.  The  usefulness  of  past 
data  depends  on  how  stationary  parameter  values  are. 

There  are  two  ways  of  using  past  measurements  to  improve  the  behavior 
of  the  controlled  system  in  the  future:  (a)  the  indirect  approach,  and 
(b)  the  direct  approach.  In  the  indirect  approach,  the  past  data  is  used 
to  estimate  the  parameters  of  the  system  explicitly  (see,  e.g.  [l],  [2], 
[3]).  In  the  direct  approach,  past  measurements  are  used  in  the  feedback 
control  law  without  going  through  the  intermediate  parameter  identifying 
step  (see,  e.g.  [■*],  [3],  [6]).  When  the  indirect  approach  is  used,  the 
control  laws  may  stimulate  the  system  to  improve  future  identification 
accuracy  (active  or  dual  control)  or  may  be  passive  where  learning  is 
accidental  (passive  control) . 

The  dual  control  law  (indirect/active)  is  the  optimal  solution  for 
controlling  a  system  with  partially  known  parameters.  However,  its 
computational  complexity  prevents  real-time  implementation  even  for  very 
3imple  systems.  The  open-loop  feedback  optimal  control  law  (indirect/ 
passive)  also  requires  a  large  computation  time  and  is  difficult  to 
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implement  :or  real  systems.  Because  ji  che  complexity  jf  tnese  icr.enes, 
global  stability  is  often  difficult  to  analyze.  If  the  parameters  mange 
quickly,  their  estimates  based  on  past  data  may  be  obsolete  by  tne  time 
enough  daca  is  collected  to  ensure  reasonable  estimation  accuracy.  In 
addition,  it  is  not  fruitful  to  spend  control  energy  in  improving  the 
accuracy  of  estimates  of  changing  parameters. 

The  direct  approaches  do  not  estimate  che  parameters  of  che  system 
explicitly,  instead  the  past  measurements  are  used  directly  in  a  feedback 
control  loop.  If  controllers  based  on  a  direct  approach  have  fixed 
structure,  parameter  insensitive  (rather  than  adaptive,'  control  lavs  resu 
The  performance  of  a  parameter  insensitive  controller  does  not  improve 
with  time  even  if  the  parameters  are  constant,  unlike  an  adaptive  con¬ 
troller  where  average  performance  improves  as  more  accurate  estimates  of 
parameters  are  obtained.  Usually,  however,  parameter  insensitive 
controllers  are  simpler  to  implement  and  analyze. 

This  paper  deals  with  techniques  for  designing  parameter  insensitive 
control  laws.  Previous  work  in  this  area  was  done  by  Cruz,  et  al  [■*], 
Hadass  and  Powell  [o],  and  more  recently  by  Kleinman  [b]  und  Padilla 
and  Cruz  [7].  All  these  techniques  either  modify  the  gains  on  the 
state  variables  resulting  from  the  quadratic  performance  index  or  include 
quadratic  functions  of  sensitivities  of  state  variables  to  parameter 
variations  in  the  performance  index.  The  sensitivities  of  the  per¬ 
formance  measure  to  variations  in  parameters  is  usually  decreased  by 
increasing  the  feedback  gains,  thus  speeding  the  dynamics.  These 
methods  do  not  use  past  measurements  and  may  therefore  only  produce 
only  a  marginal  reduction  in  parameter  sensitivity  in  many  cases. 
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Our  approach 


for  designing  parameter  insensitive  controliers  is  sigmfi- 
canciv  different  from  the  previous  approacnes.  '.v'e  first  show  that  control 
laws  based  on  linear  combinations  of  current  state  vector  are,  in  general, 
nonootimai  with  unknown  parameters.  The  control  laws  are  then  extended 
i.a)  to  include  linear  combinations  of  past  states  and  inputs  in  addition 
to  linear  funccion  of  current  state  and  (b)  to  include  control  gains  which 
are  periodic  functions  of  time.  Both  of  chese  controllers  provide  much 
better  svscam  performance  chan  modern  controllers.  It  is  shown  by  using 
simple  first— order  linear  svstem  examples  t.nat  systems  for  wnich  stability 
cannot  oe  assured  over  tne  range  of  parameter  variations  using  conventional 
controllers  (.based  on  using  linear  combinations  of  states  for  control  inputs; 
can  be  stabilized  by  using  each  of  the  controller  structures  proposed  by  us. 
Another  important  element  of  our  parameter  insensitive  controllers  is  that 
tney  may  be  used  for  large  as  well  as  small  variations  in  parameters. 

Similar  results  are  derived  for  parameter  insensitive  filters  aisc. 

This  paper  is  organized  as  follows.  Section  1  describes  two  criteria 
of  optimality  for  parameter  insensitive  controllers.  Section  3  shows  that 
the  constanc  gain  controllers  are  ncnoptimal  within  the  class  of  linear 
controllers  when  parameters  are  not  known  exact!1-'.  Sections  -  anc  5  discuss 
two  forms  of  parameter  insensitive  filters,  whi.e  Section  o  extends  the 
concepts  to  a  parameter  insensitive  filter.  Section  7  presents  some  important 
extensions,  followed  by  an  example  in  Section  8  and  summary  and  conclusions 
in  Section  9. 

2.0  CRITERIA  OF  OPTIMALITY 

«ec  us  suppose  that  the  control  .aw  is  to  be  cesigr.ee  to  minimize  3 
oerrormance  index  which  is  i  function  of  output  mu  ccntroi  time  mist  ones. 
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* 


iiace 


he  value  of  che  performance  Index  depends  on  unknown  parameter 


j  for  a  controller  with  fixed  structure,  rcr  large  parameter  variations 

i 

we  may  define  several  criteria  of  optimality.  Two  of  these  are: 

*  i,l)  Worse  Case  Design:  la  this  case,  che  control  law  is  selected 

to  provide  the  minimum  cost  for  Che  worst  value  of  o,  i.e.,  the  control 
law  is  selected  such  chat 

max  min  J (b ,  u) 

-Z'.  u 

This  is  also  called  a  minisax  design.  The  problem  can  be  considered  to  be 
a  two-person  game  in  which  the  designer  selects  the  control  law  and  nature 
selects  the  parameter  b  from  the  feasible  sec 

(2)  3esc  Average  Design:  The  control  law  is  selected  to  minimize  the 
average  value  of  che  cost  function  over  re.'.,  i.e., 

min  <f  [j(b,u)] 
u  b?.l 

where  <?  is  the  expectation  operator. 

3.0  NONOPTIMA L 1 TY  OF  THE  CONSTANT  TAIN  CONTROLLER 
l  When  che  parameters  of  a  linear  system  are  known  exactly,  and  the 

.  performance  index  is  of  the  form  of  equations  i2)  and  <3  ,  the  control 

input  is  a  linear  function  of  the  state  variables  and  the  control  gams 
reach  steady  state  values  as  the  time  horizon  is  expanded.  In  this  section, 
we  show  that,  in  general,  linear  constant  gain  controllers  are  none pt ima 1 
from  the  class  of  linear  controllers  if  the  parameters  nav  lie  anwnere  in 
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a  specified  region.  For  me  sake  of  slap  lie  icy,  we  will  prove  this  ror  a 
first-order  continuous  system  with  two  oossiole  sets  or  parameter  values 
^extension  of  the  results  to  higher  order  system  with  many  possible  para¬ 
meter  values  is  conceptually  scraightrorward) . 

Consider  a  linear  system  wich  two  possible  descriptions 


:<1  =  a1x1  +  V  +  " 

x_  *  a?x^  +  b^u  4-  w 


(1) 


Consider  the  performance  index 


_  1  rT  ,  2  z 

^  (qx  -r  u 

“  J  a 


)  dt 


(5) 


and 


or  : 


a  control  law  of  the  form 
u  =  cx 

he  intensity  of  process  noise  is 
he  state  under  the  two  descriptions,  follow  the  equations 

h  =  2lai  T  bic)Xi  *  v 

X,  »  2  (a,  b.,c)X,  * 


and  X,  anc  X.  are  the  covaria 


6) 


and 


io  c")X  dt 


The  Hamiltonian  for  the  best  average  design  is 


1-1  2 

H  *  u  '  —  ( q  4-  c  i  X  -r  •  .  ( 2 1  a  .  +b  .  c  i  X  -  •'  • 

.  ,  2  i  i  ill 

1*1 

b.o’1  and  a,  b,c'  are  both  negative,  a  staabv  state 

followine  eauati'r.s 


r.ces 


obtained  is  the  solution  jf  the 


To  check  the  opcimality  of  the  sceady  state  solution,  we  check  the  second 
order  necessary  conditions  by  linearizing  the  state  equations  and  the 
performance  index  about  the  steady  state  solution 


"  -  2  ( a  1  mb  t  c  ) 

0 

"  2b, X,  ' 

f  G  ~ 

i.  x 

0 

-2  (a1mb_ic)_ 

.2b,X,. 

(12) 


3y  per 
vector 


H* 

xc 


forming 
,  it  is 
5Z  ■ 


<52J 


Lc_bi'i’  -:+02 ’ 


3  suttaoie  set  or 


easily  shown  that 

(F-GH_1H^  )  :Z 
cc  xc 


H 

CC 


xrx2 


transformations  on  the 
[9]  equations  (12)  may 
■  Gcc 

T  -1 

-H  H  H  0 
xc  cc  xc 

0  H 

cc 


(13) 

perturbation  scat 
be  written  as 


dt 


( 1  i ) 


To  ietaraine  if  5  J  can  be  negative,  we  use  Villens  i.9_  results  on  the 
oositivitv  of  the  algebraic  Xiccati  equation  and  follow  1  orocedure  sim.i 


to  the  one  used  by  Speyer  [lo]  for  showing  the  optimality  of  airtrart  or- 
•I'iliems  shows  that  a  necessary  anti  sufficient  condition  for  ■  to  ne 


positive  is  that 

Y  <_~)  =  H  -  G*(-j 
cc 


be  positive  for  all 

YU)  =  H  1 
cc 


T  -1  T  *  r  -1  -IT 

— F  +H  a  “O')  K  H  ri.  <  <  -F-GH  *H* 

CC  CC  AC  CC  SC  cc  xc 

For  scalar  H  ,  Y  ( „)  nav  be  simDi 
cc 

* 

_ act _ 

(H  +a)  (H  +'i") 
cc  cc 


(15) 


led  to 


(16) 


where 


a  =  H1  (j^-?r1G 

xc 


(l~) 


Clearly  YO  is  negative  for  some  -  iff 


i'H** 

cc 


>  1 


(13) 


or 


Re  (a)  < 


cc 


(19) 


where  Re (a)  is  the  real  part  of  a.  It  is  easily  shown  that 


Re  (a)  =  h 


(a. 


b.,-a,,b  )“•.  (a^+b^c)  (a^+b^c) 


+  2i(a  b,-a  b^a,!),  )  (b^-b,)  (a,+a?-!-(b,+b,,)c 


c) 


CO) 


where  n  is  a  positive  constant.  Since  (a,+b,c)  and  (a7+b7c)  are 

negative,  ReO  can  be  made  as  negative  as  desired  by  selecting  a 

.2 

suitably  large  j  as  long  as  a^b7  ^  a7b,  •  Therefore,  b"J  can  be 
made  negative  and  the  constant  gain  feedback  solution  is  nonoptimal. 

It  is  interesting  to  note  that  if  a  “a,  and  b,=b7,  ReU)=0,  and  a 

constant  gain  controller  is  optimal. 
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Theorem:  '.Cher.  the  parameters  or  a  linear  svsrem  ire  nor  known  e::i 

i  controller  using  constant  gain  feedback  on  me  current  scare  is  in 
general  r.onoptimai  in  the  class  of  linear  controllers  even  in  the  stead 
state . 

The  proof  for  a  first-order  system  has  been  given  above.  It  can 
easily  be  extended  for  multivariable  system  (examining  the  system  in 
modal  form  is  one  approach).  Ve  demonstrate  using  a  simple  example  cha 
a  periodic  controller  may  provide  an  improved  overall  performance  than 
a  constant  gain  controller. 

uxamp^e :  Periodic  control.  consaoer  a  rirsc— order  examoie 

x(k+l)  =  ax(k)  -  u(k)  +■  v(k)  (21) 

where  a  is  known  to  lie  between  0  and  2.  Clearly  the  controller  u  = 
cannot  stabilize  the  system  for  ail  values  of  a. 

Consider  a  periodic  control  law  with  period  2  and  sains  c,  anc 
c?.  This  gives 

x(k-KL)  =  (a  +  c, ) x ( k )  +•  wi.k) 

(22) 

x ( k+2 )  =  (a  +  c,)x( k-f-1 )  -  w(k-t-l) 

The  covariance  of  the  state  will  alternate  between 


11 


1 

q  [  1-*-  <  a-c  ? )  "  ] 

•» 

1 - 1 a-t-c  ) “ i a-cn ) 
1  _ 


q[l+(a+c  )"] 

x,  = - T1 - -  (23) 

l-(a+c  )  (a+c,) 

and  scab ii  icy  requires 

(a+c. ) (a+c,,)  <  L  ( 2  A ) 


Stability  can  easily  be  guaranteed,  for  example,  by  selecting  c^  =  -0.5, 
c,  =  -1.5.  If  c^  =  c?  (constant  gain  controller),  the  closed-loop  system 
cannot  be  stabilized  for  all  possible  values  of  a. 

i  i 

Suppose  the  performance  index  is  ('<_+u~)  .  The  worst  case  design 
gives  c,  =  -0.37,  c^  =  -1.83  (the  worst  values  of  a  are  0  or  2).  3etter 
design  may  be  obtained  by  using  a  periodic  control  with  a  longer  time  period. 


4.0  LINEAR  OPTIMAL  PERIODIC  CONTROLLERS 

Having  'shown  in  the  previous  section  that  constant  gain  linear  controllers 
are  nonoptimal  for  systems  with  unknown  parameters,  we  investigate  controllers 
with  periodically  varying  gains.  Consider  the  system  with  state  equations 
(1) ,  performance  index  of  equations  (2)  ana  (3)  and  a  control  law  of  the  form 


u(k)  =  c,  x(k)  ( 25 ) 

k ,  31 

where  kjm  is  k  modulus  m  and  m  is  the  period  of  the  control  law. 

The  closed-loop  system  equations  are 

x(k+l)  =  (A+BC,  )x(k)  +  w(k)  (26) 

k  |  m 

and  the  state  covariance  follows  the  equation 

X(k+1)  =  (A+BC,  .  )X (k) (A+3C  ,  )  +  W  (27) 

x|m  k|m 

X (k)  =  E(x(k)xT(k) )  (23) 

The  periodic  system  is  stable  if  and  only  if  all  eigenvalues  of 
m 

II  (A+3C.)  (29) 

1-1 

are  less  than  one.  This  condition  shows  in  a  qualitative  way  the  reason 
for  the  effectiveness  of  the  periodic  controllers.  If  a  constant  gain 
controller  is  used  the  eigenvalues  of  (A+3C)  must  be  smaller  than  zero 
to  ensure  stability.  The  selection  of  a  set  of  gain  C  ^  in  a  period 
controller  increases  the  flexibility  in  the  choice  of  the  feedback  control 
structure  Co  ensure  stability  in  particular  and  to  provide  a  more  desirable 
performance  in  general.  It  is  quite  easy  to  see  that  the  eigenvalues  of 
the  matrix  of  equation  (29)  may  be  made  smaller  than  one  for  much  larger 
variation  in  elements  of  A  and  3  than  the  eigenvalues  of  the  matrix 
(A+3C) .  This  result  agrees  with  the  theorem  proved  in  the  previous  section. 

Having  selected  the  periodicity  m  of  the  controller  the  optimal 
periodic  feedback  gains  may  be  computed  as  follows.  Even  though  a  time 
varying  periodic  gain  may  be  computed  we  will  only  consider  the  steady  state 
case.  If  the  closed-loop  system  is  stable,  i.e.  all  eigenvalues  of  the 
matrix  of  equation  (29)  are  smaller  than  one,  the  periodic  steady  state 
covariances  follow  the  equation 
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Worse  Case  Design.  The  worst  case  controller  is  obtained  by  solving 
equations  (30)  and  (34)  and  selecting  3  and  C  according  to  the  following 
expression 

(  1  m  T  ) 

max  min  '  -  I  (H+DC.)  Q(H+DC.)X.' 

c1,C2,...,C;a  l  i=l  1 

m 

+  Z  Ai{(A+B.C)Xi(A+B.C)T  +  W  -  Xi+1)  (36) 

Therefore,  we  must  seek  a  saddle  point  solution  to  solve  for  3  and  C^, 
i=l,2,...,m.  The  gains  must  be  constrained  such  that  the  closed-loop  system 
is  stab la¬ 
in  the  next  seccion  ,  we  analyze  another  new  type  of  linear  controller 
structure  which  decreases  parameter  sensitivity  of  the  closed-loop  system. 

5.0  CONTROLS  LINEAR  FUNCTIONS  OF  PRESENT  AND  PAST  STATE  AND  INPUT  VARIABLES 
When  the  parameters  of  a  system  are  known  exactly,  the  present  value  of 
the  state  vector  contains  all  the  information  about  the  history  of  the  system. 
If  the  parameters  are  not  known,  however,  the  past  measurements  may  contain 
significant  information  about  the  nature  of  the  system.  All  adaptive  schemes 
make  use  of  past  measurements  to  obtain  an  improved  control  input.  Nonlinear 
functions  of  measurements  are  used  to  estimate  parameters  in  the  indirect 
approaches  and  are  used  directly  in  other  adaptive  control  schemes.  'Then  the 
parameters  are  estimated  they  are  used  in  feedback  control  calculations. 

Even  for  linear  open-loop  system  dynamics,  these  approaches  lead  to  complex 
nonlinear  closed- loop  systems.  As  mentioned  in  Section  1,  even  the  stability 
of  such  systems  may  be  difficult  to  ascertain  under  all  possible  conditions. 

Because  past  mesaurements  contain  information  about  system  parameters, 
they  may  be  used  in  parameter  insensitive  control  schemes.  It  is  simplest 


co  use  past  measurements  directly  in  a  linear  feedback,  controller.  This 
scheme  excracts  a  large  portion  of  total  information  available  from  the 
past  measurements  and  provides  a  simple  control  law  whose  properties  arc 
easy  co  analyze  and  whose  stability  regions  may  be  easily  derived. 

We  again  consider  a  system  and  a  performance  index  of  equations  (l)-(3) 
and  a  control  input  which  is  a  linear  function  of  past  states 
m 

u(k)  =■  Z  C  x(k-i+l)  (37) 

i-1  1 

Then  the  closed-loop  system  follows  the  equations 

m 

x (k+1 )  =  (A+BC  )x(k)  +  T.  BC.x(k-i+l)  (38) 

i=2  1 


If  we  define 


A'" 


A+BC.  BC„ 

1  i 

I  0 

I 


BC3  ... 


BC 


0  ...  I  o 


(39) 


WOO 


0  0 
0  0 


(40) 


x' (k)  = 


x(k) 

x(k-l) 

x(k-m+l) 


(41) 
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X'(k)  =  E[x”  (k.)  x'T(k)] 


X'(k+1)  -  A  X'  (k)  A  +  W'  (43 

The  gains  C./s  may  be  selected  as  time  varying  functions  or  as  constants  to 
optimize  the  performance  index  of  aquation  (2).  To  keep  control  laws  simple 
we  will  discuss  the  case  where  C/s  are  constant,  nonperiodic  functions. 

In  terms  of  X' (k) ,  the  performance  index  of  equations  (2)  and  (3)  is 
written  as 


J  =■  Tr  l  (H'+D'C')  Q'(H'+D'C')X' 
k=l 

where 


K  0 


a'  - 1  : 


0  0 


,  c'  =  [cx  c2 


Q  0 


0  0 


0  0 
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The  constant  Ch's  are  based  on  steady  state  system  operation.  A  steady 
state  solution  to  equation  (44)  is  obtained  if  and  only  if  all  poles  of  the 
matrix  A",  defined  by  equation  (45)  are  within  the  unit  circle.  The  steady 
state  value  of  X*  is  then 

X'  -  AWT  +  W'  (43) 


and  the  performance  index  becomes 

J  -  Tr  (H'+D'0TQ'(H'+  D'C')X' 


Worst  Case  Design.  In  the  worst  case  design  the  performance  index 


becomes 

J  =  max  Tr(H'+D'c')TQ' (H'+D'OX' 

Se  2 

Adjoining  the  constraint  of  equation  (48)  using  Lagrange  matrix  A',  it 
is  easily  shown  that 

A'  =  a'TAV  +  (H'+D'C")Tq"(H"+D>C') 


and  the  gain  C'  is  obtained  from 

max  min  Tr  [  (H'+D'C')TQ' (H'+D'C')X'+  A"  (A"X'A'T+W'-.X')  ]  (52) 

9<i2  C  " 

Best  Average  Design.  The  performance  index  for  best  average  design  is 

J  =*  E  TrOr+D'C')T(T(H'+D'cOx'  (53) 

9 

Again  using  the  Labrange  matrix  we  must  use  equation  (51)  together  with 
E  [X'(H'+D'C')Tq'D'  +  X'A>TA'B']  =  0  (54) 


Best  Average  Design  for  Small  Parameter  Variations.  Let  us  suppose  that 
there  is  a  single  parameter  3  which  varies  in  a  small  region  about  the 
nominal  value  3q.  Then  equation  (54)  may  be  written  as 


3^r  ■> 

i,  (0  )  +  itt  A0  +  — 3  A 6~  ....  =0 

°  36  30 2 


X' (U'+D'C')Q'D'  +  X'A  A'B' 


If  che  parameter  does  not  change  too  much,  we  can  neglect  all  terms  of  order 
higher  chan  two.  Lee  V  be  the  variance  of  0.  Then  equation  (56)  is 
approximated  as 

,2, 

5(9  )  +  -2-4  V  =  0  (58) 

°  302 

32" 

— 3-  is  computed  by  differentiating  appropriate  quantities  and  is  related  to 
30Z 

the  information  matrix  of  parameters  3.  Thus  the  approach  of  Kleinman  and 
Rao  (12 ]  is  a  special  case  of  the  proposed  approach. 

Selection  of  m.  The  selection  of  m  is  based  on  the  conflicting 
requirement  of  a  simple  control  law  and  reduced  storage  requirements  against 
a  smaller  performance  index.  The  marginal  improvement  often  decreases  as 
m  is  increased.  In  practice,  it  is  best  to  start  with  a  small  m  and  then 
to  increase  it  if  the  overall  performance  is  unsatisfactory. 

Special  Case.  As  mentioned  before,  all  previous  formulations  for 
parameter  insensitive  controllers  are  based  on  using  only  a  linear  function 
of  the  current  state  in  feedback.  Our  formulation  of  this  section  may  be 
specialized  to  previous  approaches  by  removing  the  superscript  prime  from 
variables  in  each  of  the  equations.  It  is  obvious,  however,  that  use  of 
past  state  variables  will  lead  to  a  lower  performance  index. 
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6.0  LINEAR  PARAMETER  INSENSITI'/E  FILTER 


The  Ideas  presented  In  the  previous  section  to  lesion  i  parameter 
insensitive  controller  may  be  extended  to  parameter  insensitive  filters.  ,vn 
optimal  filter  for  the  system  ot  equations  (1;  and  i 3;  is  given  ov 

x(k-M;k)  =  A.x(k,k)  -r  3^u(k) 

x(k  k)  »  x(kjk-l)  +  Kv(k) 

v(k)  *  y(k)  -  H.x(kjk-l)  -  D.u(k)  (59) 

where  x(k;£)  is  the  estimate  of  x(k)  based  on  measurements  till  time 
The  gain  SC  is  selected  by  solving  a  Riccati  equation  „1 3  ...  ""he  suoscritt 
represents  the  value  of  Che  quantity  used  in  the  filter.  Clearly, 


x (k+1  j  k)  =■  Af  (I-tCH-)x(kik-l)  +  (A-AfKH)x(k) 
+  (B-AfXD)u(k)  +  w(k)  -  Kv(k) 

whe  re 

x(k|k-l)  -  x(k)  -  x(kik-l) 

A  =*  A-Ae 

i. 


(60) 


(61) 


and  w(k)  and  v(k)  are  processes  and  measurement  noise  sources  with  co- 
variances  V  and  V,  respectively.  If  A,  3,  H,  and  D  are  known 
accurately,  second  and  third  terms  on  the  right-hand  side  of  equation  (6) 
drop  out  and  the  gain  K  computed  from  the  Riccati  equation  gives  the  optima 
state  estimate.  When  the  parameters  of  a  system  may  lie  anywhere  in  a  region 
the  estimation  error  is  likely  to  be  much  higher  than  normal  if  the  gain  K 
is  based  on  a  particular  Riccati  equation  solution.  In  fact,  the  error  co- 
variance  may  be  larger  than  the  state  covariance  (indicating  a  zero  estimate 
for  the  state  vector  is  better  than  the  filtered  estimate) . 
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Selection  of  control  laws  of  the  forms  of  equations  «3">  and  .  '2  adds 
new  poles  to  the  system  and  may  change  the  location  of  closed-loop  system 
zeroes.  The  ability  to  move  the  zeroes  of  the  closed-loop  svstera  is  partly 
responsible  for  the  increased  capability  to  minimize  effects  or  parameter 
variations. 

B.  A  Reformulation:  The  control  law  of  equation  (72;  be  re¬ 

form-  ated  as 

*  * 

uOc)  =  _  C  .a  (k-i-rl)  I.  *3  ' 

i-1  1 

where 

★  *  *  *  * 

x  (.k+1 )  =  A  x  (k)  +  S.u(k)  +  K_x(k) 


leading  to  the  transfer  function 


where 


u 


(I 


:  c .  2 

L=  1  1 


.1  -1  *  -  T  *  -1 

(Z  *I-A  )  -3  ) 


Z  ^x(k)  =  x(k+l). 

Equation  (71)  leads  to  tne  following  transfer  function 


-i 

u  3 

:  I  g.z1 

c  +  z  c.z1 

.  i  i 

O  _ ,  1 

L  1=1 

*  _1 

-A  i  *K*x  ',75) 


Clearly  the  two  formulations  may  be  made  equivalent. 

*  *  * 

A. ,  B_,  and  are  selected  in  a  suitable  manner  to  optimize  the 

★ 

performance  index.  Note  x  (k)  may  be  of  lower,  higher  or  the  same  order  as 
the  state  vector.  The  approach  proposed  by  Padilla  and  Cruz  [*]  is  a 
special  case  of  the  general  formulation  presented  here.  Padilla  and  Cruz 
try  to  design  optimal  filter/controller  combinations  to  minimize  sensitivity 
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Examo le :  We  cake  Che  same  system  and  show  c:ie  use  of  oast  feedback 

to  achieve  insensitivity 

x<k-’-l)  -  ax  (it)  +  bulk;  -i-  w(k)  1  ~  5 

Case  1:  b  known  exactly  and  assume  it  equals  one.  Let 

■a(k)  =  c ,  x  ( k )  +  c-x(k-i)  C'6) 

i  - 

the  region  of  stability  for  various  choices  of  c1  and  : 0  is  shown  in 
Fig.  1.  For  nonzero  c7  and  requirement  of  closed-loop  stability,  the 
maximum  allowable  variation  in  parameter  a  doubles  as  compared  to  zero 
o,  (conventional  controller).  Also  shown  in  the  figure  are  curves  of  constan 
mean  square  state  excursion,  X.  If  V,  =1  and  W.,=0,  and  parameter  a  varies  with 

L 

a  +  la,  oocimal  c  =-a  and  optimal  c^  is  shown  in  Fig.  2  for  best 
o  —  1  o  2 

average  and  worst  case  designs.  The  sensitivity  of  the  response  may  be 
further  reduced  by  incorporating  x(k-2)  in  control  u(k) . 

Case  2:  b  not  known  exactly.  This  controller  may  be  used  to  decrease 
sensitivicy  to  variation  in  b  also.  Fig.  3  shows  maximum  allowable  vari¬ 
ation  in  control  gain  b  as  a  function  of  a  and  la.  The  effect  of 

o 

not  using  x(k-l)  in  the  feedback  is  to  reduce  allowable  variation  in  b 
(see  Fig.  3). 

9.0  SUMMARY  AND  CONCLUSIONS 

This  paper  discusses  passive  controller  structures  which  minimize 
sensitivities  of  closed-loop  system  to  variations  in  system  parameters. 

These  controllers  do  not  possess  any  learning  properties,  therefore,  their 
performance  does  not  improve  with  time.  However,  they  are  simple  to  implement 
and  their  stability  regions  may  be  determined  using  well-known  techniques. 

We  developed  two  linear  parameter  insensitive  controllers,  one  based 
on  a  periodic  function  of  state  variables  and  the  other  on  a  linear 
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combination  of  past  states  and  inputs. 


'Jsing  past  state  ana  inputs  in  the 


control  lav  is  consistent  with  adaptive  controllers,  which  are  based  on 
nonlinear  functions  of  past  variables.  Leads  and  lags  used  in  classical 
control  laws  are  also  special  cases  of  the  parameter  insensitive  controllers 
suggested  by  us  in  the  paper. 

The  concepts  used  in  formulating  parameter  insensitive  controllers  are 
extended  to  the  design  of  parameter  insensitive  filters.  When  the  parameters 
3re  known  exactly  and  an  optimal  filter  is  used,  the  innovations  are  white 
[ll],  but  with  unknown  parameters  the  innovations  are  in  general  correlated. 
Past  innovations  are  therefore  used  in  addition  to  the  current  innovations 
to  update  the  state  estimate  at  each  measurement.  A  periodically  varying 
filter  gain  matrix  may  also  be  used  to  minimize  the  efface  of  parameter 
variat ions . 

In  feedback  control  problems  with  noisy  measurements,  a  parameter 
insensitive  filter  may  be  combined  in  series  with  a  parameter  insensitive 
controller  to  give  a  minimally  sensitive  compensator.  The  performance  index 
for  the  overall  system  may  sometimes  be  optimized  directly  to  give  a 
compensator  using  past  as  well  as  current  measurements  and  inputs  for  feedback. 

Parameter  insensitive  controllers  and  filters  presented  in  this  paper 
extend  the  structure  of  controllers  for  completely  known  systems  to  partially 
known  system.  logical  manner.  Our  formulation  provides  a  natural  bridge 

between  classical  an-  modern  controllers  as  well.  Classical  controllers 
based  on  incorporating  appropriate  leads  and  lags  are  often  significantly 
less  sensitive  than  modern  controllers,  but  do  not  possess  the  same  optimality 
properties.  All  classical  controllers  are  subsets  of  our  linear  parameter 


M, 
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insensitive  controllers  and  optimality  properties  can  easily  be  imposed 


using  our  approach.  This  fact  may  explain  some  of  the  well  known  properties 
of  classical  controllers. 

The  advantages  of  our  controllers  are:  (a)  they  are  simple  to  design  and 
implement,  (b)  a  desired  tradeoff  nay  be  made  between  optimality  and  complexity, 

(c)  the  stability  properties  and  stability  regions  are  easily  determined, 

(d)  the  structure  of  the  controller  does  not  change  as  long  as  the  parameters 
vary  within  the  prescribed  set,  and  (e)  the  controller  works  well  when 
parameters  change  relatively  fast.  A  disadvantage  is  that  there  is  no 
learning,  which  may  be  important  when  parameters  remain  constant  over  long 
periods.  This  is  often  a  minor  restriction,  however,  because  parameter 
insensitive  controllers  are  required  most  when  parameters  do  in  fact  change 


with  time. 
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FIGURES 


Region  of  Stability  and  Curves  of  Constant  X  for  a  Controller  3ased 
on  Linear  Combination  of  Past  States 

Control  Gain  c,  as  a  Function  of  Variation  in  a 

Allowable  variation  in  concrol  gain  b  as  a  function  of  a  and 

=  o 

variation  Isq  in  a  for  which  closed-loop  stability  may  be  guaranteed 
(broken  line  denoces  feedback  on  current  state  only  and  solid  line 
denoces  insensitive  controller  using  feedback  on  one  past  state). 
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BEST  AVERAGE 
DESIGN 


APPENDIX  B:  FREQUENCY-SHAPED  COST  FUNCTIONAL^:  AN 
EXTENSION  OF  LQG  DESIGN  METHODS 

The  1 i near -quadra t ic - Gauss i an  (LQG)  method  for  feedback 
control  design  is  extended  to  include  frequency  shaped  weighting 
matrices  in  the  quadratic  cost  f unc tional.  This  extension 
provides  a  means  to  meet  classical  design  requirements  with 
automated  computational  procedures  of  modern  control  theory.  A 
design  algorithm  to  optimize  frequency  shaped  cost  functionals 
requires  definition  of  new  states  and  the  solution  of  a  modified 
LQG  problem.  Four  examples  are  presented  to  demonstrate  frequency 
shaping  methodology:  (1)  aircraft  in  lateral  wind,  (2)  an 
industrial  crane,  (3)  vibration  control  in  helicopters,  and  (4)  a 
controller  based  on  a  reduced  order  dynamic  model. 


I.  INTRODUCTION 


Modern  control  theory  method  of  1 inear-quadrat ic - gaus s ian 
(LQG)  design  is  based  on  linear  model,  Gaussian  noise  and  quadratic 
performance  index.  The  state  variable  model  is 
x  =  Fx  +  Gu  +  w,  x(0)  =  xQ 

0  <_  t  <  tf 

where  x  is  nxl  state  vector,  u  is  qxl  control  vector  and  w  is 
nxl  Gaussian  noise  with  mean  zero  and  intensity  Q.  The  control 
u  is  selected  to  minimize: 

t  r 

d  =  fQ  (xtAx  +  uTBu)dt  +  xT(t£)Sfx(t£)  (2) 

A  is  a  positive  semi- definite  matrix  and  B  is  a  positive  definite 
matrix.  The  optimal  control  law  is  linear1. 

u*(t)  =  C ( t ) x ( t )  (3) 

C(t)  is  obtained  explicitly  by  solving  the  following  Riccati 
equation 

S  =  FtS  -  SF  -  A  +  SGB'1GTS,  S(t£)  =  Sf  (q) 

C(t)  =  -B*1GTS  (5) 

This  equation  must  be  solved  backward  in  time.  The  control  law 
is  dependent  on  system  definition  matrices  F  and  G  and 
performance  index  weighting  matrices  A,  B  and  S£ .  The 
control  gain  C(t)  is  a  function  of  time. 
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As  I  *  *»  ,  the  solution  to  Eq  .  -1  r  e  a  c h  *.*  •*  a  >  t  eacv  state 
’.'he  steady  -  state  solution  satisfies  the  a  1  get  ra : c  ?  i  c  c  a  t  ;  e-.ua- 
t i on  ARE 

-FTS  -  SF  -  A  ♦  SriB'1GTS  =  0 

Ihe  feedback  control  signal  requires  all  s  v  *  t em  states 
which,  ir.  general,  are  not  measured.  The  states  nay  be  estimate 
hv  using  an  observer  with  measurements  y  t  .  which  are  ccntan.na 
ted  with  noise  v 

y  =  H  x  ♦  v  "  j 

x  =  F  x  ♦  G  u  K  v  v  -  H  x  ,  'Sj 

In  equation  ,  x  is  an  estimate  of  the  state  vector  and  K  is 
selected  to  minimise  estimation  error  based  or.  the  process  and 
the  measurement  noise.  The  control  law  is  modified  to  use  the 
estimated  state  vector  instead  of  the  real  state. 

u  =  Cx  (9) 

It  is  clear  from  the  above  that  the  modern  control  theory 
is  based  on  optimisation  methods  with  3  specified  model  form. 

The  feedback  system  behaves  well  under  toe  following  conditions. 

i 1  The  model  is  valid  for  all  values  of  inputs  and 

states.  In  addition,  the  dynamics  is  well  described 
at  all  input  and  state  frequencies. 


(2)  The  filter  design  also  assumes  that  the  dynamics 
is  known  equally  accurately  at  all  frequencies 
(The  filter  uses  values  of  F,  G  and  H  matrices 
explicitly  in  addition  to  the  gain  K  which  also 
depends  on  the  state  definition  matrices).  This 
may  make  the  combination  of  the  filter  and  the 
control  law  extremely  sensitive  to  errors,  see  also 
Dov le“ . 

(5)  The  optimality  of  the  filter  is  strongly  dependent 
on  the  accuracy  of  noise  statistics. 

Classical  control  design  methods  are  based  on  frequency 
domain  descriptions  of  systems  and  often  account  for  model  un¬ 
certainty  at  high  frequencies  by  the  use  of  gain  and  phase 
margins  that  are  larger  than  those  needed  if  the  model  were 
certain.  The  system  dynamics  may  be  described  in  terms  of  the 
output-input  transfer  function  (s  is  the  Laplace  variable) 

y(s)  =  T (si  u ( s )  (10) 

Bode  plots,  Nyquist  diagrams  or  other  frequency  domain  methods 
are  used  to  develop  a  feedback  compensator.  In  general,  the 
control  law  is  of  the  form 

u(s)  =  Cc(s)y(s)  (Hi 

For  single- input ,  single -output  systems,  the  compensator 
Cc(s)  is  obtained  through  the  use  of  root-locus,  Bode  or  Nyquis 
plots.  The  solution  procedure  is  complex  for  multi-input,  mult 
output  systems.  Rosenbrock^,  McFarlane4,  and  others  have 
done  much  work  in  this  area,  though  the  practical  application 
of  these  techniques  is  still  quite  difficult. 


II.  PROPERTIES  OF  LQG  AND  CLASSICAL  CONTROLLERS 


The  LQG  controllers  have  a  fixed  structure  described 
previously.  The  closed  loop  transfer  function  is  of  the  form 
(,~  is  the  frequency) 

x  ( j  ui)  =  (jul  -  F  -  GC)  w(ju)  (12) 

y(ju))  *  H(ju>I  -  F  -  GC)  ^  w  ( j  u, ) 

It  has  been  shown  by  Anderson  and  Moore 3  and  Athens  and 
Safonov  that  the  controller  has  a  o0°  phase  margin  and  50* 
to  infinite  gain  margin.  The  phase  and  gain  margin  properties 
provide  for  a  constant  phase  error  in  all  channels  and  for 
individual  variations  in  gain.  The  system  may  be  extremely 
sensitive  if  two  gains  change  in  opposite  directions.  In  addi¬ 
tion,  the  phase  and  gain  margin  property  does  not  relate  directl 
to  parameter  sensitivity.  When  a  filter  is  used  for  state  esti¬ 
mation,  the  gain  and  phrase  margin  properties  are  no  longer 
valid.  In  addition,  the  filter  dynamics,  dictated  by  specifie 
noise  characteristics,  may  be  too  fast,  leading  to  interaction 
with  unmodeled  terms.  A  very  definite  attention  must  be  given 
to  the  problem  of  filter  design  when  estimated  states  are  used 
m  LQG  controllers.  The  transfer  functions  of  the  closed  loop 
system  when  a  filter  is  used  for  state  estimation  is  given  by 

x(jw)  =  [juI-F-GC(juI-F-GC*KH)  *  1k'H  ]  '  1w  ( i  u)  (13) 

y(jw)  =  Hx(ju) 

Note  that  F+GC-KH  may  have  poles  in  the  right  half-plane,  lead¬ 
ing  to  right  half-plane  closed-loop  zeros,  which  sometimes  cause 
additional  sensitivity. 

The  closed-loop  transfer  function  has  n  poles  and,  at 
most  n-1  zeros.  Therefor0,  the  frequency  response  may  drop  off 
as  slowly  as  1  at  high  frequencies  where  the  model  is  not  valid 


(higher  frequency  modes  are  often  neglected) .  This  kind  of 
response  may  be  undesirable  in  many  closed- loop  systems.  In 
addition,  the  controller  has  a  tendancv  to  push  all  responses 
to  higher  frequencies  resulting  in  more  stringent  actuator  and 
sensor  constraints.  These  LQG  properties  are  a  consequence  of 
the  optimization  problem  which  assumes  the  model  is  valid  at 
all  frequencies  and  for  all  deviations  of  the  state  variable. 

The  classical  controllers  when  properly  designed,  work  well. 
The  general  form  of  the  compensator  matrix  Cc(s)  provides 
enough  flexibility  for  a  variety  of  control  design  requirements. 
The  major  problem  is  the  difficulty  in  obtaining  the  control 
law.  The  closed  loop  transfer  function  between  the  output  and 
the  noise  is 

yCj*)  =  [I  +  T(jw)Cc(j.)]'1TK(j.)w(iu) 

where  T  Cj~)  is  the  noise  to  output  transfer  function. 

LQG  methods  with  frequency- shaped  cost  functionals  are 
developed  to  incorporate  good  classical  control  properties  in 
modern  control  design  methods,  such  that,  automated  computation¬ 
al  procedures  may  be  used  to  compute  control  laws  which  now 
require  more  difficult  graphical  approaches. 

The  next  two  sections  present  the  theoretical  background 
and  a  design  algorithm  for  the  frequency- shaping  approach. 

III.  FREQUENCY- SHAPING  OF  COST  FUNCTIONALS 

To  understand  the  concept  of  frequency  shaping,  it  is 
necessary  to  write  the  standard  LQG  cost  functional  of  Eq.  (2) 
in  the  frequency  domain.  With  infinite  time  horizon  and  no 
weighting  on  the  final  state,  the  cost  functional  may  be  written 
in  the  frequency  domain  using  Parseval's  theorem 
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J  =  1/2  /“]x*(ju))Ax(j(iO  +  u*(ju)Bu(jw)  [dui  (15) 

where  *  implies  complex  conjugate.  Clearly,  in  this  formulation, 
the  weighting  matrices  are  not  functions  of  frequency,  i.e.,  the 
state  and  control  excursions  at  all  frequencies  are  considered 
equally  unacceptable.  In  many  systems,  on  the  other  hand,  in¬ 
puts  in  the  neighborhood  of  a  particular  frequency  are  not  desir¬ 
able  because  of  poor  sensor,  actuator,  or  model  characteristics 
at  that  frequency.  Historically,  this  constraint  of  constant 
weighting  at  all  frequencies  has  resulted  because  of  the  diffi¬ 
culty  of  shaping  the  weighting  functional  with  frequency  in  the 
conventional  time  domain  LQG  formulation.  Representation  of  the 
cost  functional  in  frequency  domain  provides  a  clue  to  the  use 
of  frequency  shaping  ideas  in  modern  control  theory  techniques. 
Matrices  A  and  B  in  Eq.  (15)  may  be  made  functions  of  frequency 

to  give  a  generalized  cost  functional 

J  =  h  /"[x*(joj)A(ju))x(jw)  +  u*  (ju)B(jw)  u(jw)]  cL  (16) 

A(jw)  and  B(ju)  are  Hermitian  matrices  at  all  frequencies.  It 
appears  that  a  solution  can  be  guaranteed  if  B(jw)  is  positive 
definite  and  A(jw)  is  positive  semi -def inite  at  all  but  a  finite 
number  of  discrete  frequencies  (though  this  is  not  a  necessary 
condition) .  It  should  be  pointed  out  that  even  under  these  con¬ 
straints,  the  solution  may  not  be  easy  to  find  and  in  fact  may  not 
even  be  causal.  The  total  class  of  weighting  functions  for  which 

a  causal  solution  may  be  found  will  be  subjects  of  future  research. 
Polynomial  weighting  matrices  A  and  B  are  studied  in  the 

following  section. 


“>  "> 


IV.  CONTROL  LAW  DESIGN 


If  the  weighting  functions  A(jw)  and  B(j<u)  are  assumed 

to  be  rational  functions  of  squared  frequency,  a  systematic 

control  design  procedure,  may  be  developed  for  positive  semi- 

definite  A(juj)  and  positive  definite  B(jw).  This  is  not  a 

serious  limitation  because  a  wide  variety  of  functional  forms  may 

be  approximated  by  ratios  of  polynominals .  To  develop  a  specific 

control  design  procedure,  it  is  further  assumed  that  A(ju)  has 
rank  p,  and  B(jiu)  is  positive  definite  with  full  rank,  q 


and 


A(jw)  =  P1(jco)P1(ju)  (l-) 

B(ju>)  =  P*(ju))P2(ju;)  (IS) 

P7  are  pxn  and  qxq  rational  matrices.  Define 
P1(jw)x  =  x1  (19) 

P2(jw)u  =  u1  (20) 

If  P  (jw)  is  a  ratio  of  polynominals  in  jw  and  the 
number  of  zeros  does  not  exceed  the  number  of  poles,  eq. 

(19)  may  be  written  as  a  system  of  differential  equations  with 
output  x  . 

zl  =  F121  +  Glx 


Hlzl  +  Dlx 


Cl) 


D^  is  zero  if  the  number  of  poles  is  at  least  one  more  than 
the  number  of  zeros.  Equation  (20)  mav  also  be  written  in  terms 

of  a  differential  equation,  again  if  the  number  of  zeros  does  not 
exceed  the  number  of  poles. 

-  -i  =  F  7  z  7  +  G~,u 


1 


u  =  H7z7  +  D?u 


(221 


The  dynamic  Eq.  (1)  and  the  cost  functional  (16)  may  now  be 
written  in  terms  of  an  extended  state  vector. 
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(24) 


Defining  appropriate  vectors  and  matrices,  Eqs.  (23)  and  (21) 
become 


X  =  F:X  ♦  GXu 


J 


xV 


The  control  law  is  obtained  by  solving 
algebraic  Riccati  equation 


(25) 

(26) 


the  following  modified 


S  F- 


F1  S 


A  +  (S  G 1+  X)  B'1  (SG1+  X)T  =  0 


and 

u  =  B" 1  (S  G1  +  N)  TX  (28) 


Equation  (28)  is  written  in  an  equivalent  form  as 


u  =  CrK  +  C2z1  +  C.t2 


(29) 


74 


The  generalized  r^mroller  i:ruct|i!  e  is  then  shown  in  Figure  1. 
This  controller  has  the  form  of  a  dynamic  compensator.  The 
transfer  function  between  the  control  input  u  and  the  state 


(jw)  =  (jul  -  F1)  1G1x 

~2  =  (jul  *  1G?u  i 

[I  -  C-  (jul  -  F2)  _1G2]  uCju.)-  tc2  +  C,  ( j  u>I  -  F  )  "1G1]  xCi-; 


uCjw)  -  [I  -  C,( jul  -  F2)  '1G2J'1[C1  +  C,(ju  I  -  Fx)  -1G1]  x(ju) 

The  compensator  may  be  expressed  in  an  equivalent  form 

shown  in  Figure  2.  Note  that  the  transfer  functions  between  x^ 

and  u  and  u^  and  u  are 

x1(ju)  -  [h:  (jul  -  F1)'1G1  *  D1-u  (ju)  (33) 

Ul(ju>)  -  ;H2  (jul  -  F2 )  ‘  1GZ  +  D,"u  (ju)  (34) 

Therefore,  the  poles  of  P,  (ju)  and  P,(ju)  show  up  as  compensator 
pcles  and  zeros,  respectively  and  directly  influence  the  closed 
loop  transfer  function.  This  provides  a  direct  relationship  between 
the  form  of  the  frequency  dependence  shaping  used  in  the  cost  func¬ 
tion  and  the  structure  of  the  overall  compensator. 

If  the  number  of  zeros  in  P^(ju)  exceeds  the  number  of  poles, 

Eq .  (21)  must  be  modified.  For  example,  if  P.(ju)  has  one  more 
zero  than  poles,  Eq .  (21)  may  be  written  as 


:1  =  F1Z1  +  G1X 


x  =  z  ^  +  D^x  +  D^x 
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Using  Eq .  (1),  this  may  be  written  as 


:1  =  Fl=l 

+  G^  x 

x  1  »  H ,  z . 

+  D,x  +  D, 

1  1 

1  1 

’  Vl 

+(d1+d11f) 

( Fx+Gu) 
+  DllGl 


C  3  6  J 


Therefore  when  C  j  oo)  may  have  more  :eroes  than  poles,  a  more 
general  form  for  Eq.  (21)  is 


z  =  FlZl  ♦  Grx 

kl-i 

x1  =  H1z1  +  D.x  +  Z  D.  1  u1"1^ 

1  L  1  i=0  1 


(37) 


Where  is  the  ith  derivative  ofu.  Equations  (  22  ),  (2  3),  and 

(24)  are  also  modified  similarly. 


V.  EXAMPLES 


Example  1:  Control  of  aircraft  during  landing  in  constant 
lateral  winds.  The  equations  of  an  aircraft  are 

x  =  Fx  +  Gu  +  Pw  (38) 

T  T 

x  =  [v,p,r,d>  ],  u  [<5a,5r]  and  w  is  constant  lateral  wind 
(v,p,r  and  4>  are  lateral  speed,  roll  rate,  yaw  rate  and  roll  angle; 
oa  and  5r  are  lateral  inputs).  Control  law's  based  on  standard 
LQG  methods  are  of  the  form 

u  =  Cx  (39) 

and  for  any  finite  C  give  a  constant  steady  state  error. 
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Constant  biases  correspond  to  errors  at  :ero  frequency.  To 
reduce  the  steady  state  error  in  v  to  :ero,  an  additional  term 
of  the  following  forms  should  be  added  to  the  performance  index 


(40) 


Def ining 

.  _  v 


k 


V 


(41) 


the  additional  term  in  the  performance  index  is  .  This  is 
integral  control  I 

There  is  a  similar  problem  in  many  systems  where  the  rates 
of  change  of  control  inputs  affect  system  dynamics  leading 
to  control  saturation.  An  example  is  the  control  of  spacecraft 
attitude  using  control  moment  gyros.  To  avoid  saturation, 
additional  terms  of  the  following  form  may  be  added  to  the  perfor¬ 
mance  index 

K  (■*: ) 

uj  ” 


This  is  input  integral  control  to  avoid  saturation. 

Example  2:  Design  of  a  control  system  for  an  overhead  crane 
(Figure  5)  The  equations  of  motion  can  be  approximated  for  small 


•  a  (x,-x,) 

mx,  =  mgsinS  %  mg  2  1 ; 

£ 

mg 

n^x,  -  f  -  mgsin9  ~  f  -  — j—  (x7  -  x^ 


(43) 

(44) 
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Consider  a  case  in  which  m  =  =  1,000  lbm,  Z  =  32  ft  and 

g  =  32  ft  sec  “ .  The  open  loop  eigenvalues  are  at  0.0,  0.0,  *v/2j 
The  problem  is  to  design  a  regulator  to  control  the  position  of 
mass  hanging  from  the  crane  without  producing  serious  residual 
oscillat ions . 

A  standard  cost  functional  is  of  the  form 
1  T  i  i  2 

J  1  Lira  y  /  (a^x“+a-,xt  +  u“)dt  ( J  5 ) 

T  —■»  0 


A  good  response  is  obtained  for  a^  =  1  and  a.,  =  0.  The  time 
histories  of  x^ ,  x.,  and  u  are  shown  in  Figure  4.  Notice  that 
the  hanging  mass  continues  to  oscillate  at  about  Vz~  rad  sec 
Improved  response  can  be  obtained  by  placing  some  penalty  cn 
velocities . 

The  most  improvement  was  obtained  by  using  frequency  shaping 
methods.  Since  there  is  significant  oscillation  around  v^rad 
sec  ^ ,  it  is  useful  to  include  a  term  of  the  following  form  in 
the  performance  index. 


i  i 


This  term  is  included  with  a  unity  weighting.  The  resulting 
controller  transfer  function  is 


u(s)  = 


.64s  +  1.1  + 


.85  +  .38s 


2.2s  -  2.1  + 


s  “  +  2 

.85  .  3  S  s 

- - n - 

s“  +  2 


Xx(s) 


x0  (s) 


(47) 


Time  history  of  responses  starting  from  an  initial  condition 
error  of  lm  are  shown  in  Figure  5.  Notice  that  the  maximum 
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control  input  is  the  same.  Major  differences  in  the  applied 


control  input  occur  throughout  the  trajectory.  At  the  end  of 
10  sec.,  the  amplitude  of  oscillation  is  an  order  of  magnitude 
smaller  than  it  is  in  Figure  4.  This  shows  the  effectiveness  of 
frequency  shaping  methods. 

Example  3:  Vibration  control  of  a  helicopter.  Helicopters 
suffer  from  significant  vibration  at  NC ,  2NT.  ...  etc.,  discrete 
frequencies  where  N  is  the  number  of  blades  in  the  rotor  and  .1 
is  the  ratational  frequency .  In  the  exampl e  considered  here,  the 
equations  of  motion  describing  the  behavior  of  a  helicopter  near 
the  vibration  frequency  are  assumed  to  be 
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(48) 


w  vertical  speed  (ft/sec) 

30  coning  angle  of  the  blades  vrad) 

5coll  i5  the  collective  pitch  control  inchc 
- y i b  is  thc  oscillating  disturbance  .force) 
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1 


I 


consider  vibration  at  20  rad  sec 

5yir  =  °-2  cos  (20  r) 
where  t  is  time  in  seconds. 


i 


(*9) 


The  way  to  avoid  vibration  is  to  minimite  response  at 
u  ,  =  20  rad  sec  ^ .  A  frequency  shaped  cost  functional  of  the 
following  form  may  be  selected. 

J  -  r  ”  Qz  i p  Cjo>)  !  2  !w  (jui)l  2  ♦  R  :coU(j^)'  2  d.  (50 

where 

PC:  uj)  =  — (5i) 

4  00  -  u)  +  j  4 0 


And  ■  is  a  design  parameter.  The  cost  functional  is  converted 
into  a  time  domain  formulation  by  adding  the  following  states. 


This  equation,  combined  with  the  helicopter  model,  specifies  the 
plant  description  to  be  used  with  the  performance  index  describee 
in  equation  (4).  The  augmented  state  description  is: 


l 

i 


r 


0 . 0339 


0  .  1 1  S  o 


I  i  i 


The  optimal  control  law  is  of  the  form: 


coll 


=  CW 


K 


i  - 


An  implementation  of  this  vibration  controller  is  shown  in  Figure 
6.  The  vibration  level  at  a  point  in  the  Helicopter  after  the 

controller  is  switched  on  is  shown  in  Figure  ~  u=.01).  Even 
though  ^  can  be  set  to  cero,  a  small  value  of  1  W3S  used  to 
minimize  sensitivity  to  parameter  variations. 

Example  4:  Control  cf  Structures.  Dynamics  of  structures 
is  described  by  a  set  of  partial  differential  equations.  Finite 
order  model  representations  are  valid  only  at  low  frequencies. 

We  show  how  frequency  chaping  methods  are  used  to  design  control 
laws  for  models  valid  only  at  low  frequencies. 

Consider  a  free  standing  pyramid  where  we  wish  to  control 
line-of-sight  errors  at  the  apex  (Fig.  S) .  The  size  of  various 
elements  are  given  in  Tables  1  and  1.  We  have  an  eight  mode 
model  representing  the  behavior  of  the  system  below  1  H„.  These 
actuators  and  three  sensors  are  available  to  actively  increase 
damping  ratios  of  specific  modes  which  affect  open  line-of-sight 
errors.  All  modes  have  3n  open  loop  damping  ratios  cf  C.?‘ . 


*  This  example  was  developed  by  Dr.  Robert  Strunce  and  his  assoc¬ 
iates  at  Charles  Star!-:  Draper  Laboratorv,  Cambridge,  '-'ass. 
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Since  the  model  is  valid  only  at 


desirable  to  minimize  control  and  state  activities  above  1  Hz. 
Ihis  may  be  achieved  by  a  control  weighting  matrix  which  in¬ 
creases  with  frequency,  for  example: 

B  (J  )  =  Ab(d"  +  lj02)  (35) 


diagonal  matrix  with  elements  b^,  b,,  and  b - •  By  selecting 


the  control  weighting  at  1.5  H. 


five  times  higher 


at  0.5  Hz.  The  control  weighting  is  written  equivalently  as: 


luTV)  (. 


3  1  . .  n  1 


U  +  w  U  =  V 


The  1 ine-of- sight  control  requires  increased  damping  in  modes  1,5, 
i  and  3.  This  is  achieved  bv  state  weighting  of  the  form 


0.5  ( q  i 


*  V4  *  q,  +  a. 


*  ta.3  +  q5'  US’ 

•a.  is  the  modal  deflection  correspond ing  to  the  i  th  mode. 

b,  ,  b,  and  b.  are  ail  set  at  .01. 
a 

Table  5  shows  the  closed  loop  eigenvalues  for  the  S  mode  desig 
model  and  the  15  mode  evaluation  model.  Because  of  the  low  control 
activity  at  high  frequencv,  the  damping  ratios  of  unmodeled,  high 
frequency  modes  have  small  perturb3t ions . 

When  a  comparable  design  was  carried  out  without  frequency 
shaping  of  state  and  control  weighting  matrices,  the  perturba¬ 
tions  in  the  unmodeled  eigenvalues  was  so  large  that  some  of  the 
modes  were  unstable. 


Discuss  ion :  We  have  shown  some  simple  applications  of 
frequency  shaping  methods.  The  method  is  much  more  powerful 
and  is  being  applied  to  many  complex  problems. 

A  reviewer  provides  an  interesting  interpretation  of  the 
controllers  of  Examples  2,  3  and  4.  Example  2  and  3  give 
"virtual  tuned  vibration  absorber"  while  Example  4  gives  a 
"virtual  actuator  dynamics". 


VI.  CONCLUSIONS 

Frequency  shaping  of  quadratic  cost  functionals  can  pro¬ 
vide  additional  flexibility  in  developing  control  laws  using 
modern  control  theory  techniques.  For  good  design,  the  part¬ 
icular  selection  of  weighting  functions  depends  on  fa)  model 
validity,  fb)  actuator /sensor  charact er i s t ic s ,  (c)  control 
objectives,  3nd  (d)  disturbance  spectrum.  Examples  were  presente 
to  demonstrate  the  effectiveness  of  the  proposed  anproach. 

The  design  procedure  for  f requencv  shaped  weighting  on  state 
and  control  law  is  implemented  by  linear  elements  with  memory. 

The  frequency  shaping  concept  has  been  extended  for  state 
estimation  with  models  V3 1 i d  within  a  finite  frequency  range  and 
with  finite  bandwidth  sensors.  These  ideas  may  be  also  be  used 
in  parameter  identification. 
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Table  1 

Cross  Sectional  Areas 


TRUSS  ELEMENT  NO. 

! 

AREA 

! 

1 

j 

100 

|  2 

100 

!  3 

100 

!  4 

100  { 

;  5 

100 

!  6 

100  ! 

7 

1000 

S 

100  0 

9 

1000 

10 

1000 

11 

100 

i  12 

1 _ 

100 

Table  Z 


Structural  Node  Coordinates 


!  NODE 

X 

Y 

-  ! 

I  1 

0.0 

0.0 

10.165 

1  ? 

j  w 

-5.0 

-2.837 

2.0 

!  3 

5.0 

-2. S87 

2  .  0 

4 

0.0 

5.7735 

2.0 

|  5 

-6.0 

-1.1547 

0.0 

j  6 

-4.0 

-4.6188 

0.0 

n 

4.0 

-4.6188 

0  . 0 

\  3 

6.0 

-1.1547 

0 . 0 

1  9 

2.0 

5.7755 

0.0 

i  r\ 

1  O 

L 

-2.0 

5 .7735 

0.0 
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APPENDIX  C: 


THh  ACCOMUDAT I  UN  UP  LARGE  PLANT  PAkA.MbTfck 
VARIATlUNs,  BY  sTAl'E  EEED8ACK 


I'h  1  s  section  examines  the  stability  properties  of  a 
cLosea-loop  system  under  a  state  feedback  control  law  ana  in  the 
presence  of  parameter  variations  in  tne  system's  matrix.  It  is 
shown  that  if  the  variations  are  of  a  special  class,  the  stability 
properties  of  the  closed-loop  system  will  De  insensitive  to  these 
variations.  This  class  of  systems  can  be  mace  even  larger  if  a 
certain  given  matrix  is  staDle.  This  analysis  distinguishes 
systems  where  linear  state  feedback  control  laws  wor<  adequately 
and  systems  where  more  detailed  adaptive  control  laws  are  requirea. 


Consider  the  following  linear,  continuous-  t  ime  uyr.  ami 
system  described  by 

x  ( t ;  =  A  x  1. 1 -*•  3u  ;  t ; 


where  x 

n 

-:K‘ 

,  uejR"1  are  the 

s  ta 

te  and 

inpu 

t  veotc 

r  <  r  e  $  a  e 

t  i  v  e  i  y  , 

A 

0 

is  an  n<.n  mat 

r  ix 

and  3 

:  s 

an  n  ■  m 

m *  n 

matrix  o 

-  r 

uil  rank.  Furtne 

rmo  r 

s  it 

is  as 

sumed  t 

h at  t h e 

matrix 

remains  constant 

but 

the  paramet 

e r s  of 

the  ti  a  t  r 

are  sub j 

ca  r 

to  changes  over 

a  p 

o  s  .s  i  b  1 ' 

v  i  a  r 

ge  kno 

wn  ran^e 

These  variations  of  parameters  are  either  induced  ur.aer  :p 

ing  conditions  or  are  simply  cue  to  error  introduced  iurin 

modeling  process.  In  either  case,  the  parameters  of  the  m 

A  are  bounded  within  the  known  range  of  oossibie  oarame 
o  • 

variations,  namely 

Armin;SA  I  A  max 

where  the  inequalities  are  taken  component - w i s e  .  Another 

wav  of  representing  the  matrix  A  is  bv 

3  o 

A  =  A  ”  A , 

0  1 


whe  i  e 

A  is  the  design  (or  nominal 

matrix  and 

A .  1st h 

X 

mat  x 

of  deviations  from  the  design 

values.  It 

s  ii o u  *. u  be 

po  inte 

d  out,  however,  that  while  A 

is  taken  as 

a  known  , 

s  tant 

matrix,  A,  can  chance  over  t 

he  complete 

r  a  n  g  .■  e  r  p 

s  ib  1  e 

parameter  variations  and  its  v 

alue  will  no 

t  r*.  c  t*  i  a  *  *  ■" 

known 

to  the  des i gne r . 

T 
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matrix  as  in  Eqs.  (2.2)  and  (2.5).  .  .n  one  find  a  constant-gain 

state  feedback  control  law  that  will  keep  the  closed- loop  system 
stable  over  the  complete  range  of  possible  parameter  variations0 
For  the  problem  to  be  solved  successfully,  it  should  specify  the 
conditions  under  which  it  is  solvable  and,  also,  identify  the 
class  of  matrices  for  which  the  stability  properties  of  the 

closed-loop  system  will  be  insensitive  to  parameter  variations. 
Ideally,  one  would  like  to  design  a  closed-loop  system  that  will 
be  insensitive  even  to  large  parameter  variations. 

The  next  section  will  present  a  particular  method  for  the 
derivation  of  the  state  feedback  control  law.  The  structure 
offered  by  this  method  will  prove  very  useful  in  answering  the 
questions  posed  above. 
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III.  STATE  FEEDBACK 


A  new  way  toward  the  output  feedback  problem  was  developed 
in  [9],  introduced  in  [10],  and  later  extended  to  state  feedback 
problems  in  [11].  This  new  approach  was  termed  the  moda lator 
approach  since  the  solution  of  the  pole-placement  problem  is 
obtained  via  a  particular  dynamic  system  that  provides  informa¬ 
tion  on  the  modal  structure  of  the  closed-loop  system  under  the 
particular  state  (or  output)  feedback  control  law. 

For  convenience  and  easy  reference,  the  main  modalator 
result  for  pole-placement  by  state  feedback  will  be  given  below, 

[11]. 

Note :  Prime  will  denote  transposition,  i.e.,  A " 

is  the  transpose  of  A  . 

THEOREM  3.1  (MAIN  MODALATOR  RESULT) 

Consider  the  linear,  controllable,  continuous  -  time*  dynamic 
system  described  by 

x(t)  =  Ax(t)  +  Bu(t)  ,  xaRn,  ue]Rm  ,  m<n  .  l3.1) 

Then,  one  can  always  find  a  stable  (n-m)x(n-m)  matrix  F  with 
arbitrary  eigenvalues  and  an  (n-m)*m  matrix  H  such  that  the 
(n-m)xn  transformation  matrix  Q  that  solves 

j 

QA'  -  F'Q  =  HB'  (3.2)  J 


*This  procedure  is  applicable  to  di screte - time  dynamic  svstems 
as  well. 
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allows 


m  n-m 

Then,  the  following  will  hold: 

U)  F'  a  QA'E  (3.4) 

H  ■  QA'D  .  (3.5) 

(ii)  If  the  state  feedback  control  law  is  given  by 
u(t)  =  Kx(t)  where 

Ka-D'i.A+aI),a>0  (3.6) 

then  the  closed- loop  matrix,  A  ,  given  by 

A  ■  Q'E'A  -  iBD'  (3. ') 

will  have  n-m  of  its  eigenvalues  equal  to  those  of 
the  matrix  F  and  the  remaining  m  closed- loop 
eigenvalues  will  be  clustered  at  (-a)  ,  i.e., 

t  (A)  d  [  \  I  a  L  \  -^  ct ) m  :F-\In.m;  =  0.  (3.3) 


Proof :  The  statements  regarding  the  solution  of  Eq .  (3.3)  , 

the  inversion  of  Eq.  (5.3),  and  the  expressions  given  in  Eqs. 
i5.4  and  (5.3)  are  well  known,  (see,  e.g.,  [12]  and  [13]). 

If  one  is  using  the  state  feedback  gain  matrix  given  by 
Eq.  (5.6),  then  the  closed-loop  matrix,  A  ,  is  given  by 

A  s  A*BK  "  (I-BD')A-aBD' 
but,  using  Eq.  (5.5),  this  becomes 

A  »  ■)  '  E  '  A  -  -a  3  D  "  . 
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The  closed- loop  eigenvalues  are  found  from 


ff(A)  =  |Q'E'A-aBD'-Xln|  =  |  Q'E '  (A+aIn)  -  ( \*<x)  I  |  =*  0 


This  can  be  simplified  [14]  to  yield 


tt  (A)  »  (\  +  a)m|  E'AQ'-XI 


(\+a)m  | F- A I 


n-m  ‘ 


n-m ' 


This  procedure  for  the  derivation  of  the  state  feedback 
control  law  results  in  a  special  property  that  can  be  observed 
by  examining  the  proof  of  theorem  3.1  and  will  be  stated  here 
explicitly . 

Corollary  3.1  (Separation  Property):  The  choice  of  the 
scalar  a  in  the  derivation  of  the  state  feedback  gain  matrix 
in  Eq.  (3.6)  will  not  affect  the  choice  made  for  the  n-m  eigen¬ 
values  of  the  matrix  F  in  Eq.  (3.2). 

The  importance  of  this  corollary  lies  with  the  fact  that 
once  F  is  determined,  changing  the  scalar  a  will  have  no 
effect  on  the  eigenvalues  of  F  but  will  only  move  the  m 
clustered  eigenvalues  along  the  real  axis.  This  properlv  is 
realizable  because  the  modalator  approach,  summarized  in  theorem 
3.1,  factors  the  n-dimensional  state  space  into  two  invariant 
subspaces . 

The  procedure  presented  here  involved  a  certain  matrix 
equation  whose  solution  provides  the  state  feedback  gain  matrix 
ir  a  closed-form  expression.  This  matrix  equation  is  usually 
associated  with  the  Luenberger  observer  [15],  and  it  is  of 
interest  to  examine  the  relation,  if  any,  between  the  observer 
and  the  modalator  approach.  This  relation  will  be  examined 
through  the  concept  of  duality. 


The  system  in  Eq.  (3.1)  given  by 


jx(t)  =  Ax ( t )  +  Bu ( t) 

|v(t)  =  x(t) 

has  a  dual  system  given  by 

x * ( t )  *  A'x*(t)  +  u*(t) 
y*(t)  =  B'x*(t)  . 

It  can  be  easily  verified  that  if 

u*(t)  =  - (A'+oI  )DB'x*(t) 

then  the  vector  defined  by 
eft)  A  QX*(t) 

evolves  dynamically  according  to 
:(t)  -  F'e (t)  . 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


Furthermore,  it  can  be  shown  that  the  transformation  matrix,  Q  , 
in  addition  to  solving  QA'-F'Q  =  HB'  ,  also  solves 

QA'-  F'Q  =  0  .  (3.13) 

Equation  (3.13)  is  particularly  interesting  because  it 
reveals  that  the  matrix  Q  is  a  matrix  expressible  by  eigen¬ 
vectors  of  A  that  correspond  to  the  common  eigenvalues  of  A 
and  F  .  The  modal  structure  revealed  through  Q  motivates 
the  terminology  chosen  for  the  modalator.  Furthermore,  this 
structure  explains  the  absence  of  an  output  term  in  Eq.  (3.12). 

We  will  not  elaborate  this  point  any  further;  the  interested 
reader  is  referred  to  [10]  and  [11]. 


IV.  STABILITY  PROPERTIES  AND  PARAMETER  VARIATIONS 


The  objective  of  this  paper,  as  stated  in  Section  II,  is 
the  derivation  of  a  state  feedback  control  law  that  ensures  the 
stability  of  the  closed-loop  system  even  in  the  presence  of 
parameter  variations.  Therefore,  one  would  like  to  design  an 
appropriate  state  feedback  control  law  such  that  the  actual 
closed-loop  matrix  given  by 

Aq  -  A  ♦  Al  (4.1) 

where 

A  I  A  *  BK  (4.2) 

will  remain  stable  over  all  possible  variations  in  A^  . 
Obviously,  this  problem  may  not  be  solvable  for  every  case  and, 
therefore,  one  should  indicate  the  class  of  systems  for  which 
the  problem  is  solvable. 

The  nxm  input  matrix  B  was  assumed  to  be  of  full  rank 
equal  to  m  .  therefore,  without  loss  of  generality,  we  can 
consider  the  matrix  B  as  given  in  the  form 


B 


(4.5) 


where  I  is  the  mxm  identitv  matrix.  When  B  is  not  in  the 
m 

form  of  Eq.  (4.3),  one  can  employ  a  simple  transformation  that 
yields  an  equivalent  system  whose  input  matrix  is  in  the  form  of 
Eq.  (4.3).  The  details  of  this  transformation  can  be  found  in 
Appendix  A. 


Since  the  matrix 
columns  of  B  and  Q' 


V 

;q\ 

can 


of  Eq.  (3.3)  is  nonsingular,  the 
serve  as  a  basis  for  the  n-dimens iona 1 


i 


can  be  written  as 


state  space,  hence  the  matrix  A, 

Al  -  BR  +  Q'T 

where 

R  is  an  mxn  matrix 
T  is  an  (n-m)xn  matrix. 


(.  4 . 4  J 


The  matrices  R  and  T  are  determined  by  the  (given)  range  of 
possible  parameter  variations.  The  actual  values  of  these 
matrices,  however,  are  not  known  to  the  designer. 

Using  Eqs.  (4.1),  (4.4)  and  (3.7),  the  actual  closed-loop 
matrix  is  given  by 

A  =  Q'E'A  -  aBD  '  +  BR  +  Q'T  .  (4.3) 


Since  BD'  *  Q'E'  =  In  ,  Eq.  (4.3)  becomes 

A  =  Q'E  ' (A+al  )  -  xl  +  BR  +  Q'T 
o  x  v  n  n  x 


which  can  be  written  compactly  as 


Ao  = 


b:q' 


__i 

[_E  (A+aIn"l  +  T  J 


al 


(4.6) 


A  question  of  interest  will  be  the  following  one.  Can  one 
specify  a  class  of  R  and  T  matrices  that  allows  the  stabili¬ 
zation  of  Aq  over  all  parameter  variations?  A  sufficient 
condition  for  an  affirmative  answer  is  given  below. 


Theorem  4.1:  If  the  procedure  of  theorem  3.1  is  used  in 
the  derivation  of  the  state  feedback  gain  matrix,  and  if 


T'cN(Q) 

R'eN(Q) 

where  N(Q)  denotes  the  null  space  of  Q, 
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(4.7) 

(4.8) 


i 


then  the  actual  closed-loop  matrix,  A  ,  will  have  n-m  of  its 

eigenvalues  equal  to  those  of  F  and  m  eigenvalues  equal  to 

those  of  the  matrix  (RB-al  )  ,  i.e., 

m 

ni)  =  !  (RB-al)  -  XI  j  -  I  F- A I  ;  =  0. 
o'  1  v  m'  m 1  1  n-m’ 


Proof:  From  Eq.  (4.6)  we  find 


t(a0) 


|_E'(A+aIn)+T  j  1_ 


(X  +  a) In  I  =  0. 


Using  Eqs.  (5.5)  and  (5.4),  this  becomes 


f  RB  I  RQ"  1  . 

ff(*V  =  |E"(A+aIn)B+TB|  F+aIn_m^TQ'J  *  (X+ot)In  |  =  0  1 

i  ! 

However,  if  R  and  T  are  such  that  Eqs.  (4.7)  and  (4.8)  hold, 
then 

*(A0)  ■  |(RB-c,)-XImMF-Un.J  -  0  . 

Note  that  if  Eqs.  (4.7)  and  (4.8)  hold,  then  the  eigenvalues  of 
F  will  not  change.  Equation  (4.9)  also  indicates  what  is  the 
role  of  the  positive  scalar  a  . 

Corollary  4.1:  If  the  positive  scalar,  a  ,  is  chosen  such 

that 

51  >  WRB> 


then  the  actual  closed- loop  system  will  remain  stable  for  all 
parameter  variations  satisfying  Eqs.  (4.7)  and  (4.3). 


The  results  presented  above  show  that  n-m  actual  closed- 
loop  eigenvalues  will  be  insensitive  to  changes  in  and  that 

the  remaining  m  closed-loop  eigenvalues  can  be  ensured  to 
remain  stable  by  a  proper  choice  of  a,  provided  the  R  and  T 
matrices  are  in  the  null  space  of  the  transformation  matrix  Q  . 
It  may  be  argued,  however,  that  since  one  starts  the  design  with 
a  given  range  of  possible  parameter  variations,  it  may  not  be 
possible  to  find  an  appropriate  transformation  matrix  Q  that 
will  satisfy  Eqs.  (4.7)  and  (4.3).  This  point  is,  of  course, 
valid;  and  we  will  therefore  turn  our  attention  to  identify  a 
class  of  systems  whose  given  variation  will  satisfy  Eqs.  (4.7) 
and  (4.3)  and,  therefore,  will  be  accommodated  by  a  constant- 
gain  state  feedback  control  law. 

It  is  well  known  (see,  e.g.,  [15])  that  if  one  starts  with  an 
invertible  matrix  in  Eq.  (3.3),  then  the  solution  of  Eq.  (3.2) 
amounts  to  a  pole-placement  problem  for  the  matrix  F.  Therefore, 
if  one  partitions  A  according  to 


and  i  i 


then,  if  B 


A  = 

“All 

;  An" 

}  m 

(4.11) 

_  A21 

;  A22  _ 

}  n-m 

m 

n-m 

Q  = 

Ql^n 

-m_  * 

(4.12) 

B'  = 

[V!0 

j  ,  it  can  be  shown  that 

F  =  E 'AQ  ' 

=  At  £  *  *4t  2 

(4.13) 
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and  the  purpose  of  is  to  achieve  the  desired  pole-placement 

for  F  .  Then 


B'  “ 

-1 

_ 

r  Im 

0 

- 1 

^m  !  I* 

Q 

L  Qi 

^n-m  _ 

_  -Qi  i  !n- 

(4.14) 

iVe  are  now  in  a  position  to  identify  the  class  of  matrices 

that  satisfy  Eqs.  (4.7)  and  (4.8). 


Let 


Di!r:! 


m  n-m 


[Tj  iL-] 

m  n-m 


(4.15) 


(4.16) 


then,  Eqs.  (4.7)  and  (4.8)  requires  that 
R2  =  -RjQ^ 

t2  -  -W  • 


C-i.i") 


The  matrix  can  now  be  written  as 


A:  *  BR  +  Q'T 


- 1 

po 

+ 

1  O 

i  *-* 

1  \ 

;  h 
»— * 

-'WW) 

Dmi 

L  Ti 

-TiV  J 

T1  . 
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and  using  Eq.  i4.14),  this  becomes 


A 


1 


Ti  ! 

A  a 


14.13) 


This  shows  that  as  long  as  A-,  has  the  form  given  in  Eq. 
(4.13),  the  actual  closed-loop  system  will  have  n-m  eigen¬ 
values  that  are  insensitive  to  the  changes  in  A,  and  m 

3  1 

eigenvalues  are  those  of  the  matrix 

F.  =  RB  -  a  I 
1  m 

which  in  this  case  reduces  to 

F.  =  R.  -  al_  .  '4.1?; 

11m 

The  class  of  systems  satisfying  Eq.  (4.13)  is  obviously 
quite  restricted.  However,  close  examination  of  Eq.  (4.  IS) 
indicates  a  wider  class  of  system  whose  stability  properties 
will  be  insensitive  to  parameter  variations. 

Theorem  4.2:  If  the  fn-m)x(n-m)  matrix  A,-,  in  Eq . 
(4.11)  is  stable,  and  one  selects  the  matrix  F  such  that 

f  =  a22  (4. 20) 

then  the  actual  closed-loop  system  will  be  completely  insensi¬ 
tive  to  variations  in  A-,^  ,  and  any  variations  in  .4^  can  be 
accommodated  by  choosing  a  such  that 

F 1  =  R  l 


is  stable. 


Ill 


A,,  is  stable  and  its  eigenvalues 
=  0  and ,  from  Eq .  :  4  .  1  3  ; 


where  R^  and  T ^  are  arbitrary.  As  was  shown  before,  n-m 
eigenvalues  will  be  those  of  F  and  they  will  be  insensitive  to 
parameter  variations,  while  the  remaining  m  eigenvalues  are 
those  of  F^  that  can  always  be  made  stable  by  choosing  a 
large  enough. 

Remark :  If  the  input  matrix,  3  ,  is  not  in  the  form 

given  by  Eq .  (.4.5)  ,  then  the  stability  proper¬ 

ties  of  cne  closed-loop  matrix  will  not  be 
determined  by  the  matrix  A^2  ,  but  a  differ¬ 
ent  matrix.  The  details  can  be  found  in  the 
appendix . 

At  this  juncture,  a  comment  is  in  order  regarding  the 
choice  of  the  design  matrix  A  in  Eq.  (4.5).  It  is  obvious 
that  different  choices  will  result  in  different  A^ .  If  one 
chooses  the  design  matrix,  A,  such  that 

A  =  A(min)  (4.22) 

then 

Ax  >  0  (4.23) 

where  the  inequalities  are  taken  component -wise ,  Since  the 
matrix  A,  is  non-negative,  bounds  can  be  found  on  its  eigen- 

X 

values  by  either  rows,  or  columns  sums.  In  that  case,  a  can 
be  chosen  to  exceed  these  bounds  and  then  one  bypasses  the  need 
to  evaluate  the  largest  eigenvalue  of  the  matrix  R^  ,  as  is 
required  by  Eq.  (4.10).  It  should  be  pointed  out,  however,  that 


Proof:  Notice  th: 


l  r 


are  choker,  cor  F  then 


__1_' 

T.  '  0 

f  i 
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whatever  choice  is  made  tor  the  design  matrix  A  ,  it  should  be 
such  that  will  be  a  controllable  pair. 

The  design  procedure  presented  thus  far  exhibits  two  prop¬ 
erties  that  should  be  pointed  out  explicitly.  Both  have  to  do 
with  the  dimension  of  the  input  vector.  First,  it  is  seen  that 
the  larger  this  dimension  is,  more  parameters  can  vary  without 
affecting  the  stability  properties,  and  more  variations  can  be 
taken  care  of  by  a  proper  choice  of  the  single  scalar  a  . 
Second,  the  larger  this  input  dimension  is,  more  closed-loop 
eigenvalues  will  be  clustered  on  the  real  axis.  .An  interesting 
case  occurs  when  the  state  and  the  input  vectors  are  of  equal 
dimensions,  i.e.,  m=n  .  In  this  case,  the  modalator  approach, 
as  given  by  theorem  (3.1),  seems  to  break  down  since  there  is  no 
matrix  equation  like  (3.2)  to  be  solved.  This,  however,  is  not 
the  case,  since  when  m=n  ,  all  closed-loop  eigenvalues  are 
clustered  at  one  point  on  the  real  axis.  The  state  feedback 
gain  matrix  is  then  given  simply  by 

K  =  -B'L(A+aIn)  ,  x>0  .  (4.24) 

This  results  in 

A  =  A+BK  =  - xl  (4.25) 

n 

and 


It  is  obvious  that  in  this  case  (m=n)  ,  one  can  always  find  a 
scalar  a  that  will  accommodate  any  variations  introduced 
through  A^  . 

The  conditions  imposed  by  Eqs.  (4.7)  and  (4.3)  resulted  in 
the  invariance  of  the  eigenvalues  of  F  to  changes  in  A^  . 
Then,  if  A-,-,  was  stable,  one  could  design  a  constant-gain 
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state  feedback  control  law  that  will  yield  closed-  1  jop  star:!:' 
properties  that  are  insensitive  to  parameter  var: at;  ;ns .  "".e 

class  of  systems  to  be  treated  can  be  enlarged  bv  dropping  c:r. 
dition  id. 7).  In  this  case,  in-mj  eigenvalues  of  the  closed 


.a.:' 


loop  system  will  be  found 

from 

T(Fo)d|F+TQ'-\In 

-mi  "  0 

Let 

,  _  j  V1’1)  - 

A0U,2) 

0  LA0(2,1)  ; 

A0l2,2) 

14.23, 


then,  if  there  exists  a  set  of  parameters  that  results  in 
A  (2,2)  being  stable,  then  one  can  set  F  3  A0i2,2)  ,  where 
A  (2,2)  is  the  one  particular  stable  matrix  out  of  the  pos¬ 
sible  variations  for  Aq(2,2)  .  Then,  as  was  shown  before, 
the  transformation  matrix  Q  is  given  by 

«  '  ("hn-J 

and  if 

Fq  -  F  +  T7  Id. 50) 


remains  stable  for  all  possible  T,  matrices,  then  the  closed- 
loop  system  will  remain  insensitive  to  parameter  variations. 

The  matrix  A^  will  be  given  in  this  case  by 


A 


1 


.  T  I 

L  ' 


(4.51) 


The  results  thus  far  were  developed  for  continuous  time  sys¬ 
tems;  Appendix  B  presents  the  extension  to  discrete-time  dynamic 
systems.  An  example  will  illustrate  the  suggested  procedure. 
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Example :  Consider  the  linear  dynamic  system  described  by 


x(t)  =  Aq x ( t )  +  Bu(t) 

where 


and 


B  = 


-5  <  a  <  2 


1  <  b  <  4 


Let 


5 

i\ 

/da 

0  \ 

0  <  Aa  < 

1 

J 

Ai  = 

L 

o) 

0  <  Ab  < 

Since  A^  =  -1,  and  (A,B)  is  a  controllable  pair,  we  can 
design  a  state  feedback  control  law  that  will  accommodate  all 
the  variations  in  the  parameters  a  and  b  . 


Since  a>Aa  is  required  for  the  stability  of  the  actual 
closed-loop  system,  as  established  in  Eq.  (4.  ),  let  a=6  , 

then 


K  =  -  D ' (A+al ) 


-(1 


-1)  • 


The  actual  closed- loop  matrix  is  given  by 


A 


o 


A+BK+A1 


'-6+Aa  0 

k  1  +  Ab  -1 
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As  seen  from  the  above  example,  the  changes  in  b  do  not  affect 
the  stability  of  A  ,  while  those  of  a  can  be  taken  care  of 
by  a  proper  choice  of  the  scalar  a  . 
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CCNCLUSKNS 


This  paper  presented  the  sufficient  condition  under  which 
the  stability  properties  of  the  closed-loop  system  will  not  oe 
affected  bv  parameter  variations  in  certain  suomatrices  of  the 
svstem’s  matrix.  Specifically,  it  was  shown  that  if  the  varia¬ 
tions  in  the  parameter  satisry  a  certain  rule  or,  if  a  certain 
submatrix  is  stable,  a  state  feedback  control  law  can  be  found 
such  that  the  stability  properties  of  the  closed-loop  system 
will  not  oe  affected  bv  parameter  variations. 


117 


APPENDIX  A 


Consider  the  system  described  by 
x  ( t  J  ■  Ax  (  t )  *  B  u  ( t ) 


and  where  che  nxm  matrix  3  is  partioned  by 


3  ,  |  ;  n-m 


(A.  Ij 


(a.:) 


Since  the  matrix  3  is  assumed  to  be  of  full  rank,  one  can 
3 1  wavs  renumber  the  components  of  the  state,  if  necessarv,  so 
that  3,  is  an  rum  nonsingular  matrix. 

Let 


3  "  1  o 

b  ,  b  : 1 ;  i 

.  1  '  n-m 


(A.  3) 


t  «  Tx . t  i 


A.  i) 


then  the  equivalent  system  to  Eq .  (A.l)  will  be  given  by 

: i t ;  »  Ac ; t i  *  3u(t )  (A. 5 i 

where 


A  *  TAT 


3  -  TB 


-  I 


m 

“o 


i  A .  0 


{ A. 
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Since  Eqs.  (A. 1  and  (A. 5)  are  equivalent  systems,  one  can 
analyte  Eq.  (.A.  5)  instead  of  Eq.  (A.  1  j  . 

The  matrix  A-,-,  will  be  given  by 

At  i  =  A_j  *  Aj  .  ( A .  8  J 


If  one  solves 

Q  A'  -  F'  Q  =  H  3' 
for  Q  ,  then 

QA  -  F'Q  =  HB ' 
is  solved  too,  with 

F  -  F 
V  H  *  H 
(  Q  *  Q(T') ‘ 1 


(A.  9) 


(A. 10) 


CA. 11) 
(A. 12) 
(A. 15) 


and 


r  1  *  1 

- 

-l 

r 

B  1 

B 

....  |  =  T ' 

— 

=  T' 

D 

L  Q  1 

.Q. 

(A. 14) 
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APPENDIX  B 


The  discrete-time  dynamic  system  is  described  by 

x (k+1)  »  AQx(k)  +  Bu (k)  t8-1- 

where  xe®n,  ue®m,  and  A  and  B  are  appropriately  dimensioned 
matrices . 

The  results  reported  in  section  III  are  applicable  to  the 
system  given  by  (B.l).  The  results  of  section  IV,  however,  should 

be  slightly  modified. 

Theorem  (4.2)  applies  to  discrete-time  dynamic  systems  except 
that  now  (4.21)  will  require  that  the  scalar  a  will  be  chosen 
such  that 


which  leads  to 


X(R1)  -  1  <  a  <  A ( R1 )  +  1. 

The  condition  expressed  in  (B.3)  may  not  be  satisfied  by  a 


(B.3) 


single  value  of 


if 


has  large  variations. 


REFERENCES 


1.  Cruc,  J.  B.  and  Perkins,  W.  R. ,  "A  New  Approach  to  the 
Sensitivity  Problem  in  Multivariable  Feedback  System 
Design,"  IEEE  Trans.  Aut.  Control.,  Vol.  AC-9,  July  1964. 

2.  Salmon,  D.  M. ,  "Minimax  Controller  Design,"  IEEE  Trans . 
Aut.  Control,  Vol.  AC-15,  pp.  569-376,  August  1963 . 

3.  Horisberger,  H.  P.  and  Belanger,  P.  R. ,  "Regulators  for 
Linear,  Time  Invariant  Plants  with  Uncertain  Parameters," 
IEEE  Trans.  Aut.  Control,  Vol.  AC- 21,  pp.  '05 -"08, 

October  1976. 

4.  Patel,  R.  V.,  Toda,  M.  and  Sridhar,  B.,  "Robustness  of 
Linear  Quadratic  State  Feedback  Designs  in  the  Presence 
of  System  Uncertainty,"  IEEE  Trans.  Aut.  Control,  Vol. 


AC- 22, 

No.  6, 

pp.  945-949, 

December  1977. 

Young  , 

K.  D.  , 

Kokotovic,  P. 

V.  and  Utkin,  V 

•  I-, 

,  "A 

Singular  Perturbation  Analysis  of  High-Gain  Feedback 
Systems,"  IEEE  Trans.  Aut:  Control,  Vol.  AC- 22,  .No.  6, 
pp.  951-93F^  December  i9'“t'. 

6.  Utkin,  V.  I.,  "Variable  Structure  Systems  with  Sliding 

Modes:  A  Survev,"  IEEE  Trans.  Aut.  Control,  Vol.  AC-22, 

No.  2,  pp.  212  -  222,  April  115":. 

Cruc,  J.  B.  and  Sundarajan,  N. ,  "Sensitivity  Improvement 
of  Feedback  Systems,"  Report  R-509,  Coordinated  Sciences 
Laboratorv,  University  of  Illinois,  Urbana,  Illinois, 

May  1971.' 

3.  Harvey,  C.  A.  and  Pope,  R.  E.  ,  "Study  of  Synthesis  Tech¬ 
niques  for  Insensitive  Aircraft  Control  Systems,"  NASA 
CR-2805,  April  1977. 

9.  Arbel,  A.,  "Aggregation  and  Control  in  Large-Scale  Dynamic 
Systems,"  Ph.D.  dissertation,  department  of  Engineering- 
Economic  Systems,  Stanford  University,  October  1978. 

10.  Arbel,  A.,  and  Tse,  E. ,  "Modalators:  A  New  Approach  Toward  the 
Output  Feedback  Problem,"  submitted  for  publication  to 

IEEE  Trans,  on  Aut.  Control. 

11.  Arbel,  A.  and  Tse,  E.,  "Modalators:  Theory  and  Some 
Applications,"  submitted  for  publication  to  Automat ica . 


121 


Aoki,  M.  and  Huddle,  J.  R. ,  "Estimation  of  the  State 
Vector  of  a  Linear  Stochastic  System  with  a  Constrained 
Estimator,"  IEEE  Trans.  Aut.  Control,  Vol.  AC- 12,  pp. 
432-453,  August  196~. 

Yuksel,  Y.  0.  and  Bongiorno,  J.  J.,  "Observers  for  Linear 
Multivariable  Systems  with  Applications,"  IEEE  T r ans .  Aut . 
Control,  Vol.  AC-16,  No.  6,  pp.  605-613,  December  19? 1. 

Plotkin,  M. ,  "Matrix  Theorem  with  Applications  Related  to 
Multivariable  Control  Systems,"  IEEE  Trans.  Aut.  Control, 
Vol.  AC- 9,  pp.  120-121,  January  1964”. 

Luenberger,  D.  G. ,  "Observers  for  Multivariable  Svstems," 
IEEE  Trans.  Aut.  Control,  Vol.  AC-11,  No.  2,  pp .  190-197, 
April  1966 . 

Gantmacher,  F.  R. ,  The  Theory  of  Matrices,  Chelsea,  New 
York,  Vol.  I,  1959. 


122 


1 

I 


APPENDIX  D:  CONTROL  OF  FLEXIBLE  STRUCTURES 


Flexible  structures  such  as  an  aircraft  wing,  long  beam,  or  a 
large  space  structure  (LSS)  are  representatives  from  a  class  of 
systems  with  infinite  modes.  The  actual  dynamics  are  governed  by 
partial  differential  equations  for  which  the  lower  frequency  modes 
can  be  adequately  (in  most  cases)  modeled  by  interconnected 
oscillatory  systems.  The  difficulty  arises  in  that  the  dimension 
of  this  model  is  extremely  large  (e.g.,  thousands  of  modes),  and 
the  higher  modal  frequencies  and  modes  shapes  are  poorly  known. 
However,  the  structure  of  the  model  is  definitely  known--i.e., 
positive  symmetric  matrix  second  order  equations.  The  problem 
then  is  to  find  robust  controllers  which  are  based  on  a  lower 
order  model  of  the  flexible  structure.  Two  approaches  are 
presented . 


D.l  ROBUST  INVERSE  OPTIMAL  CONTROL  FOR  FLEXIBLE  STRUCTURES 


An  inverse  optimal  control  problem  is  Formulated  to  develop 
robust  control  laws  tor  purely  oscillatory  systems.  The  optimal 
control  solution  requires  output  feedback  with  specif iea 
constraints,  leading  to  robustness  with  respect  to  unmouelea  modes 
and  a  large  class  of  parameter  variations.  The  robustness 
properties  are  proved  directly  from  known  properties  of  control 
laws  optimizing  quadratic  performance  indices.  These  control  laws 
are  useful  for  poorly  damped  flexible  structures. 


1.  INTRODUCTION 


Ever  since  Kalman  (Ref.  1)  formulated  the  inverse  optimal 
control  problem  in  1964,  the  subject  has  been  of  significant 
interest  in  control  theory  (Refs.  2'51.  The  inverse 
optimal  control  theory  attempts  to  derive  a  class  of  performance 
indices  which  are  optimized  by  a  prespecified  control  input. 

In  its  simplest  form,  the  state  and  input  penalty  matrices  of 
a  quadratic  performance  index  are  determined  for  linear  feedback 
control  of  linear  dynamic  systems  (Refs.  1  and  6). 

Past  research  has  never  exploited  the  full  potential  of 
inverse  optimal  control  theory.  Once  the  control  law  has  been 
specified,  obtaining  the  optimizing  cost  functional  itself  has 
little  direct  value.  However,  when  combined  with  known  and 
previously  proven  properties  of  optimal  control  laws,  the  derived 
optimizing  function  aids  the  understanding  of  robustness  prop¬ 
erties  of  the  feedback  control  law.  An  optimizing  function  is 
derived  for  a  specific  class  of  systems  and  is  used  to  study 
controller  properties. 

This  paper  addresses  purely  oscillatory  systems  defined 
by  all  poles  and  zeros  on  the  imaginary  axis.  Flexible  systems 
such  as  large  space  structures  (LSS)  are  examples  of  such  systems, 
because  of  their  typically  low  damping  ratios.  Optimizing  functions 
are  derived  for  control  laws  which  increase  damping  ratios  in 
such  systems.  The  optimizing  functions  are  quadratic.  Since 
the  systems  are  linear,  the  properties  of  linear  quadratic  controllers 
are  used  to  specify: 

(1)  robust  control  structures; 

(2)  the  set  of  parameter  variations  for  which  the 
controller  is  robust;  and 

(3)  the  class  of  truncated  modes  under  which  the 
closed  loop  system  is  stable. 


126 


The  inverse  optimal  control  problem  considered  here  sim¬ 
plifies  the  process  of  obtaining  optimizing  functions  for  a 
given  control  law.  It  is  shown,  for  example,  that  the  procedure 
for  obtaining  quadratic  cost  functions  is  simplified  considerably 
if  the  dynamic  systems  are  expressed  in  modal  form.  The  paper 
places  an  emphasis  on  studying  methodologies  to  utilize  the 
optimizing  cost  functionals.  A  system  for  which  these  results 
are  directly  applicable  is  that  of  a  large  flexible  space 
s  t  rue  tu  re  ( LSS  )  . 

The  special  forms  of  dynamic  models  are  given  in  Section 
2.  Section  3  presents  the  problems  to  be  considered.  Section 
4  solves  the  inverse  optimal  control  problem  and  Section  S  presents 
some  further  results  associated  with  this  solution.  Section 
.6  discusses  the  important  issue  of  robustness  of  the  closed- 
loop  system. 


2 .  MODEL  FORM 

The  dynamics  of  continuous  flexible  structures  are  described 
by  second  order  partial  differential  equations  of  the  wave  type. 
The  partial  differential  equation  is  not  amenable  to  control 
design  and,  therefore,  approximate  dynamic  models  are  obtained 
through  the  application  of  the  finite  element  method.  The  dynamic 
behavior  of  a  flexible  system  can  be  closely  approximated  by 
the  following  undamped  structural  equations 

Mn  +  Kn  =  F  (1) 


where  n  is  a  vector  of  generalized  displacements  at  various 
points  on  the  structure.  Equation  (1)  can  be  transformed  to 
the  following: 


«•  A  q 


=  Bqu, 


qeR  ,  ueR1 


(2) 
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where  the  nxn  matrix  A 

o 

(p.d) 


is  symmetric,  and  positive-definite 

(3) 


These  equations  are  known  as  the  structural  modal  equations 
since  the  modal  structure  can  be  easily  derived  for  the  system 
described  in  Eq.  (2).  Defining  a  state  vector 

x  =  [q,q]T  (4) 


Eq.  (2)  can  be  written  in  state-space  form  as 


x(t)  =  Ax(  t )  +  Bu(  t )  ,  xc]R2n,  udRm 


where 


A  = 


0  1 

l 

i 

n 

0 

B  = 


0 

B 


(5) 


(6) 


and  Bq  is  an  nxm  matrix  of  full  rank  m  (m<n). 

The  structural  properties  of  Eqs.  (5)  and  (6)  are  suitable 
for  the  description  of  the  flexible  space  structure.  Therefore, 
Eqs.  (3),  through  (6)  define  the  modal  form  to  be  considered 
in  the  inverse  control  problem. 

The  purpose  of  control  in  flexible  systems  is  to  provide 
damping  to  the  oscillating  structure.  This  is  accomplished 
through  energy-dissipating  devices.  To  implement  these  controls, 
velocities  are  measured  at  different  locations  on  the  structure 
and  converted  to  inputs  through  actuators  (rate  feedback). 

When  the  sensors  and  actuators  are  co-located,  the  m-dimensional 
output  vector  is  obtained  through 


yf  t )  =  BTx( t ) ,  y cRm 


(7) 


where  B  is  the  same  matrix  appearing  in  Eq.  (5)-. 
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Let  us  assume  for  the  moment  that  the  positive-definite 
(p.d)  matrix  AQ  in  Eq.  (6)  is  a  diagonal  matrix  (this  can 
always  be  done  through  a  simple  transformation),  i.e. 


Aq  =  diag(a1 ,  a,,  ....  a  )  ,  ai  >  O.vi  <_  i  _<n 


(S) 


Then  it  can  be  verified  (Ref.  7)  that  the  i-jt^  element  of  the 
matrix  of  right  eigenvectors  U  of  the  2n  x 2n  matrix  A  in 
Eq.  (6)  is  given  by 

1  i  =  1,2 . n 

j  =  2i-l,  2 i 


and 


u  •  . 
i) 


x . 
J 


i  =  n+1 ,  . . . ,  2n 
j  =  2(i-n)-l,  2  (i -n) 


(9) 


where 


elsewhere 


AU  =  UA 

(10) 

A  =  diag(X1,  X 2  »  • • • » 

(11) 

(  j  /a ( i+ 1 ) / Z  1  =  1 ’ 3’  ‘ 

Xi  *  „  i  =  2,4,  . 

(  -j  yai/2 

..,  2n- 1 

.  .  ,  2n 

(12) 

i  *  . 


Then,  using  Jacobi's  formula  for  eigenvalues  perturbations,  it 
can  be  shown  (Ref.  8)  that  the  use  of  rate  feedback  where  the 
sensors  and  actuators  are  co-located  will  result  in  an  eigenvalue 
shift  to  the  left  in  the  complex  plane  (increased  damping).  The 
control  law  is  an  output  feedback  given  by 
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u(t)  =  Ky  (t)  =  KBTx(t)  f 13) 

This  form  of  direct  velocity  feedback  can  be  shown  to 
dissipate  energy  and,  therefore,  stabilize  the  system  provided 
that  energy  in  any  zero  frequency  modes  is  conserved  (Refs.  9  and  10) 
The  proof  is  based  on  using  the  energy  as  a  Lyapunov  function. 

The  observation  that  the  closed-loop  system,  under  the  condi¬ 
tions  described  above,  is  stable  motivates  the  interest  in  the 
inverse  optimal  control  problem.  The  following  factors  further 
define  closed-loop  control  requirements. 

The  dimension  of  the  model  described  in  Eq.  (5)  depends  on 
the  fineness  of  the  grid  used  in  the  application  of  the  finite 
'element  method.  Usually,  however,  the  dimension  is  quite  large 
and  some  reduction  in  order  is  required  to  facilitate  numerical 
solutions.  This  reduction  in  order  will  result  in  modes  that  are 
truncated  from  the  final  model  used  for  design.  The  effect  of 
the  control  derived  from  the  reduced-order  model  on  the  truncated 
modes  is  of  extreme  importance  (control ' spi 1  lover ,  Ref. 11) 

The  solution  of  the  inverse  optimal  control  problem,  to  be  presented 
shortly,  will  be  brought  to  bear  on  this  problem  in  Section  6. 


3.  PROBLEM  STATEMENT 


Consider  the  linear  dynamic  system  described  by 

jx(t)  =  Ax(t)  +  Bu(t),  xeR2n,  ueRm 
I  y(t)  «  B Tx ( t )  ,  vcRm 

where 


A  - 

r  o  i  i  I 

B  = 

0 

i - 

i 

> 

o 

o 

1 _ 

-  V 

(14) 


(15) 
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and  the  nxn  matrix  AQ  is  symmetric  and  pos i t i ve - de f in i te < p • d ) 
i .  e  . 


A 


o 


>  0 


(16) 


Next,  consider  the  set  of  admissible  controls,  y/,  defined  by 

A  ^  x 

?/  =  (ueR  :  u  =  Ky,  Re  [  A  (A  +  BKB  )  ]  <  0}  (17) 


where  X ( ■ )  stands  for  the  eigenvalue  of  the  argument,  and  Re ( • ) 
is  the  real  part  of  a  complex  number. 

Let  the  objective  function  to  be  minimized  be  quadratic  in 
the  state  and  control  over  an  infinite  time  period,  i.e. 

•  J(x,u)  4  (x^(t)Qx(t)  +  u^(t) Ru (t ) ]  dt  (18) 


with 


Q  1  0  ,  R  >  0 


(19) 


Problem  1:  (inverse  optimal  control  problem) 

Given  an  admissible  control,  u* ,  under  what  conditions  on 
K,  Q,  and  R  is  this  control  law  optimal?,  i.e.,  when  is 


J(x,u*)  = 


min  J(x,u) 
ue  V/ 


(20) 


Problem  2: 


To  what  extent,  if  any,  is  the  solution  to  the  inverse  optimal 
control  problem  sensitive  to  truncated  modes? 

Problem  3: 


How  are  closed-loop  stability  properties  affected  by  param¬ 
eter  variations  (or  uncertainty)  in  system  matrices? 
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4.  INVERSE  OPTIMAL  CONTROL  PROBLEM 

Part  of  the  problem  presented  in  Section  3  is  that  of  finding 
for  what  state  and  input  oenalty  matrices,  Q  and  R,  respectively, 
a  given  admissible  control  also  minimises  an  LQ  problem.  If 
such  an  optimal  control  u*  exists,  then  (Ref.  2) 

u*  =  Ky  =  -R'TBTPx  (21) 

where  P  is  the  2n  x 2n  p.d  solution  of  the  algebraic  Riccati 
equation  (ARE)  given  by 

PA  +  ATP  -  PBR*  ^  BTP  +  Q  =  0  (22) 

The  inverse  optimal  control  problem  is,  therefore,  the 
following.  Given  an  admissible  rate  feedback  control  law  find 
for  what,  if  any,  weight  matrices  Q  and  R,  the  ARE  in  Eq.  (22) 
is  solved. 

Let  the  2n  x 2n  symmetric  P  and  Q  matrices  be  partitioned 
as 


P  = 


11 


r  pi  i  j  pu] 

Q  = 

Q11  1 

pT,  !  p„ 

.  i  ‘-4-- 

0  1 

[  1 

ARE  in  Eq . 

(22) ,  we  find 

AoP22  ‘  P 

i  7B  R"  1B^'  P_,  = 
X  Z  0  O  Z  b 

0 

T 

+  A  P !  „  + 

O  O  1  L 

P1?B  R '1bI  P.. 
1 2  o  o  1 Z 

=  Qn 

-IT 

R  AB*P-,- 
O  o  2  2 

p  pT  =  Q 

12r12  ^2  2 

0 

Q22 


(23) 


(24) 

(25) 

(26) 
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Since  a  rate  feedback  control  law  is  used,  we  find  from 
F.  q  .  (211  that 

3Tp  =  0  —  P,  ,£  V(BT)  i  f  ;,-Rn  :BTc  =  0  !-  ''2') 

Therefore,  ,  can  be  expressed  as 

pi:  =  \yn-J  i28) 

where 

Ny.VfB  Tq)  (2  9) 

f  =  *T  >  0  (30) 

The_  n  xn  matrix  P^-,  as  given  in  Eq.  (28)  is  symmetric  and  from 
Eq.  f30).  p.s.d.  The  latter  condition  is  reouired  to 
satisfy  Eq .  (25)  since  (-A  )  is  a  stable  matrix. 

Lemma  1 : 

For  a  dynamic  system  under  a  rate  feedback  control  law,  the 
nxn  matrix  in  Eqs.  (24)  through  (2b)  has  to  satisfy 

PI2  ■  0  (31) 

Proof : 

This  is  verified  by  observing  Eq.  (26)  where  Q.,.,  has  to  be 
at  least  p.s.d.  ■ 

Since  P^  =  0»  the  resulting  2n  x  2n  matrix  P  that  solves 
the  ARE  will  be  a  block  diagonal  matrix.  This  motivates  the  follow¬ 
ing  form  for  the  matrix  P  to  be  checked  as  a  possible  solution 
of  the  ARE. 


Let 


P  =  a 

8(BT 

-IT  T 

B)  B  ♦  NEN 

(32) 

where 

NeV(B 

T), 

NeR2n  x  ^2n'm) 

(33) 

NTN  - 

I  , 

2n-m 

(34) 

cteR1 

+ 

E  =  E 

T  > 

0  E£R^2n’m^  X  ^2n‘m^ 

(35) 

It  is  clear  that  P 

’  in  Eq.  (32)  is  symmetric. 

To 

show 

that 

it  is  also  p.d. 

we 

note  that 

PB  = 

aB 

(36) 

PN  = 

NE 

(37) 

Therefore,  B 

and 

N  are  right  eigenvectors  of 

the 

matrix  P 

corresponding  to  positive  eigenvalues;  therefore 

,  P 

is 

p.d. 

Similar  results 

are  obtained  for  left  eigenvectors. 

Using  the 

form  for  P  as  given  in  Eq.  (32) 

the 

ARE 

can  be 

written  as 

PA  + 

T 

A  P 

2  - 1  T 

-  a  BR  B  +  Q  =  0 

(38) 

Since 

'  0  ' 

B  = 

— 

A. 

-i 

I 


i 


i 
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we  find  that  to  satisfy  Eq .  (33) 


0 

1  I 

1  n 

,\l 

1  0 

M  cRn  X  fn'm) 


(39) 


where  N  is  chosen  such  that 
o 


BfN  =  0  «•  N  e  KbT) 
00  o  o' 


T 

N  N 
o  o 


n-m 


(40) 

(41) 


Let  the  symmetric  matrix  E  in  Eq .  (32)  be  partitioned  as  follows 


1  E,  "1 

_ i..: 

T  1 

E2  1  E5 


i  n-m 


(42) 


)  n 


n-m  n 


Theorem  1: 


A  necessary  and  sufficient  condition  for  P  in  Eq.  (32)  to 
solve  the  ARE  in  Eq.  (22)  is  that 


0 


0 


I  i  -IT 

0  I  a ^ B  R  XB 
o  o 


(43) 


E  = 


'arn-m| 

0  ■ 

0  * 

aA_  - 

(44) 


The  2n  x 2n  p.d.  solution  matrix  P  is  given  by 


P  = 


aA„  I  0 

.9.  J - 

0  I  ol 


(45) 
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Proof:  (.sufficiency 


0  |  aB  (BTB  )  '  1B1 

I  OOO  0 


a  A  1 
o  I 


0  I  iN  N  1 
i  o  o 


since,  from  Eqs.  (40)  and  f 4 1 ) 


T  1-1 


T  -  1  i 
B  f  B  1  B  )  |  N 

OOO  I  0 

I 


we  get 


sA  i  0 
0  I 


0  1  si 


■'  -IT  T 

Q  =  s"BR  B1  -  PA  -  A*P 


(necessity)  from  Eq .  (58) 


Q  =  a2BR_1B  -  PA  -  A  P 


Since 


1  -IT 
0  |  s“B  R  B 
0  0 


E,  I 

j 


T  T 

e;n1 

2  o 


N  E ,  1  iB  (BTB  )  " 1 BT  *  N  E.M 
o  2  |  o  oo  o  olo 


we  find 
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Q  = 


E^N  A 
Zoo 


+  A  M  E, 
o  o  2 


T  -IT 

A  [aB1  (B  B)  B 

O  O  K  0  o 


X  EvNT] 
o  1  o 


[aB  (BTB  )'LBT  +  N  E1NTJA  -E- 
L  o  oo  o  olo  o  a 


2  -IT 
a  B  R  B  -  N .  E 


-  E  _N 

o  2  2  o 


Since  Q  is  block  diagonal,  then  using  Eq.  (46)  and  Lemma  1,  we 
find 


E,  -  0 

E ,  =  a  A 
a  o 

Once  the  solution  to  the  ARE  is  obtained,  one  can  find  the 
optimal  control  as  well.  The  case  under  consideration  in  this 
paper  is  that  where  the  sensors  and  actuators  are  co-located. 
For  this  system  we  have  the  following. 

Coroll arv  1: 


When  the  output  feedback  control  law  is  chosen  as 

u  =  Ky  (47) 

where  .<  is  a  symmetric,  mxm  negative-definite  (n.d.)  matrix 
and  the  output  is  given  through 

y  =  BTx  (48) 


then 


R  =  -aK 


(49) 


IS" 


w 


m 


is  the  m  xm  control  penalty  matrix  of  the  quadratic  objective 
function 


J  =  (xTQx  +  u TRu)  dt 


(50) 


Proof : 


The  optimal  control  that  minimizes  Eq .  (50)  is  given  by 
u*  =  -  R~ 1 B  TPx 

using  the  solution  matrix  P  from  Eq .  (45)  we  find 

"  "  -1  T 

u*  =  -aR  B  x 


hence 


K  =  -aR 


5.  FURTHER  RESULTS 

In  Section  2,  we  stated  that  the  closed-loop  system  with 
colocated  sensors  and  actuators  is  stable.  The  proofs  are  based 
on  Jacobi's  formula  for  eigenvalue  perturbations  and  energy  func¬ 
tions.  The  solution  obtained  in  the  previous  section  can  be 
used  to  prove  the  closed-loop  stability  properties  of  the  system 
under  the  control  law  derived  in  Eq .  (47). 

The  stabilization  of  a  linear  system  by  pole-placement  via 
state  feedback  requires  the  system  given  by 

x  ( t )  =  A  x  ( 1 1  +  B  u  ( t ) 

to  be  fullv  controllable  iRef.  12)  (This  requirement  can  be 
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relaxed  to  that  of  stabil izability  fRef.  13)).  This  is  an  impor¬ 
tant  structural  property  that  will  be  assumed  for  the  system  under 
cons i de  rat i on . 


and  the  observability  matrix  of  the  pair  (A,L)  bv 


a 


(A,L)  = 


T  i  r  T 
L  >  A  L 


J-,2n-  ] ,  T 


•  I 


--  (A1)-'*  ‘  L 1 

then,  using  the  expression  for  L,  we  can  verify  that 
■  A_ 


( A  ,  L )  =  a 


0 


0  '  I 


n  J 


( A ,  B )  H 


1/2 


.  i  / 1 

where  the  2nm  x 2nm  matrix  H  is  defined  bv 


1/2  A 


-1/7 

R  L/~ 

1 

1 

_  _  J 

o 

0 

0 

-1/7  1 

i 

0 

t 

t 

‘  \ 
i 

0 

,  \ 

0  ‘  -  a. 

i 

r-1/2 

A  simple  application  of  Sylvester's  inequality  yields  the  required 
result.  g 

The  results  obtained  thus  far  can  be  used  to  yield  a  physical 
interpretation  of  the  solution  to  the  inverse  optimal  control  prob¬ 
lem. 

Since  the  state  vector,  x,  is  partitioned  according  to 
T 

x  =  [q  ,q]  the  optimal  performance  is  written  as 


J*  =  min[J^  +  J,] 
u 

(54) 

where 

JI  ■  «A0*  1TB0R‘1bJ  q  dt 

1 

(55)  1 

■J  ?  =  /”  uTRu  dt 

| 

(56)  * 
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The  interpretation  suggested  above  is  that  the  cost  functional 
is  the  total  energy  of  the  system.  The  term  contributed  by  the 
state  is  the  kinetic  energy  where 

E  =  qTMq 

and  the  "mass,”  M,  is  given  bv 

M  =  a"3  R_1BT  r  5  7 ) 

o  o  ; 

Also,  the  use  of  velocity  feedback  where  the  actuators  and  sensors 
are  colocated  results  in  an  LQ  problem  where  only  kinetic  energy 
is  weighted. 

The  results  of  the  inverse  optimal  control  problem,  as  dis¬ 
cussed  thus  far,  are  summarized  in  Table  1. 


0.  ROBUSTNESS  PROPERTIES 

As  mentioned  in  Section  2,  one  of  the  problems  associated 
with  modeling  large  flexible  space  structures  is  that  of  truncated 
modes.  The  matrix  describing  the  dvnamics  of  the  flexible  system, 
given  by  F.q.  (6)  has  eigenvalues  distributed  on  the  imaginary  axis. 
The  process  of  mode  truncation  is  carried  out  by  ignoring  eigen¬ 
values  that  are  far  from  the  origin  on  the  imaginary  axis  (high- 
frequencies).  Analytically,  the  problem  of  mode  truncation  is 
as  follows.  Let  the  full  order,  finite  linear  system*  be  described  by 


*  The  r e a  1  flexible  system  is  described  through  a  second  order  PDE 
wave  equation).  In  linear  form,  the  flexible  system  is  of  in¬ 
finite  dimension.  Equation  (58),  therefore,  considers  only  a 
finite  number  of  these  modes. 
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Table  1 
Summary 


PROBLEM: 

J*  ~  min  r  (x^(t)Qxi't)  +  uT(t)Ru(t))  dt 
u  -  //  0 

s.  t. 

x  ( t )  =  Ax(f.)  +  Bu(t),  xe!R^n  ,  uePm 
u ( t )  -  KB^x(t) 

K  =  <  0  (A  given  gain  matrix) 


x  =  [q,q]T,  x(0)  =  xQ 


INPUT  PENALTY  MATRIX: 

R  =  •  a  ^ ,  a>0 

STATE  PENALTY  MATRIX: 

Q  -  f-Q-i — 9 — 

0  'a2B  R"V 

L  0  0 

SOLUTION  OF  THE  ALGEBRAIC  RICCATI  EQUATION  (ARE): 


P  = 


_'o"iTn_ 


OPTIMAL  COST: 


J*  ■  xTx0  »  «(qjA0%  * 


where  the  subscript  T  denotes  the  part  of  the  system  to  be 
truncated.  The  output  used  for  control  is  given  by 


y(t)  =  B^x(t)  «■  B^xT(t),  ye  ]Rm  (59) 

The  matrix  A,  contributes  to  the  full  order  system  frequencies 

higher  than  those  contributed  by  A  and,  therefore,  can  be 

0 

truncated  to  yield  the  reduced  order  model  described  by 

uft)  (60) 

-  -  The  first  robustness  result  can  now  be  stated. 

Theorem  5: 

The  full  order  system  described  in  Eq.  (58)  under  the  control 

u(t)  =  Ky ( t ) ,  -K  =  -KT  >  0  (61) 

and  y  given  by  Eq .  (59)  will  "emain  stable  for  every  Ay  provided 

that  (A  ,  3  ),  (Ay,  By)  are  controllable  pairs. 

Proof : 


X  (t)  = 

r°  "nl 

---4--- 

x(t)  + 

0 

-A  |  0 

L  r\  1 

.  B  . 

Using  the  control  as  derived  from  the  reduced  order  system 
results  in  the  full  order  system  assuming  the  form  given  by 


2  Ct) 


Theorems  1  and  2  indicate  that  this  system  will  be  asymptotically 
stable  provided  (Ao>Bq)  and  form  controllable  pairs. 

This  result  shows  how  the  inverse  optimal  control  solution 
can  be  put  to  use  to  answer  the  question  of  control  spillover 
(Ref.  9). 

Observing  the  proof  of  Theorem  1  we  can  arrive  at  the  second 

robustness  result.  This  result  shows  how  stability  properties  of 

the  closed-loop  system  are  affected  by  parameter  variations  in 

the  svstem  matrix  A  . 

o 

Coro  1 lary  2 : 

'  "  Consider  the  dynamic  system  described  in  Eq .  (60)  and  where 

the  n  x  n  matrix  A  can  have  uncertain  parameters  that 
belong  to  a  certin  set.  Applying  the  output  feedback  control 
of  Eq .  (61)  will  result  in  a  stable  cloosed-loop  system  for 
every  A  in  the  given  range  of  possible  parameter  variation 
provided  A  is  always  symmetric  and  p.d,  over  its  entire  range. 

Proof : 

From  theorem  1. 

l 

Remark:  In  the  case  when  Aq  is  a  diagonal  n  x  n  matrix,  the 

above  result  guarantees  stability  of  the  cloosed-loop  system 
provided  the  diagonal  of  AQ  remains  positive  over  all  possible 
variations . 
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Theorems  1  and  2  indicate  that  this  system  will  be  asymptotically 
stable  provided  C  Aq  ,  Bq  )  and  form  controllable  pairs. 

This  result  shows  how  the  inverse  optimal  control  solution 
can  be  put  to  use  to  answer  the  question  of  control  spillover 
(Ref.  9J . 

Observing  the  proof  of  Theorem  1  we  can  arrive  at  the  second 

robustness  result.  This  result  shows  how  stability  properties  of 

the  closed-loop  system  are  affected  by  parameter  variations  in 

the  svstem  matrix  A  . 

o 

Corollary  2: 

'  '  Consider  the  dynamic  system  described  in  Eq .  (60)  and  where 

the  n  x  n  matrix  Aq  can  have  uncertain  parameters  that 
belong  to  a  certin  set.  Applying  the  output  feedback  control 
of  Eq.  (61)  will  result  in  a  stable  cloosed-loop  system  for 
every  AQ  in  the  given  range  of  possible  parameter  variation 
provided  Aq  is  always  symmetric  and  p.d.  over  its  entire  range. 

Proof : 

From  theorem  1. 

■ 

Remark :  In  the  case  when  AQ  is  a  diagonal  n  x  n  matrix,  the 

above  result  guarantees  stability  of  the  cloosed-loop  system 
provided  the  diagonal  of  A  remains  positive  over  all  possible 
variations . 
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D . 2  STABILITY  OF  LQG  MODAL  CONTROL  FOR  FLEXIBLE  STRUCTURES 


The  staoiiity  of  LOG  (.linear  observer  based  feedoaCK.)  moaai 
control  of  lirge  flexiole  structures  is  examined.  It  is  shown 
that  the  effect  of  the  ooserver  and  the  spillover  (residual  modal 
interaction  with  primary  modes)  is  eauivalent  to  a  finite  gain, 
linear  function  perturbing  an  LQ  (linear  primary  mooal  state 
feedback)  control.  Recent  results  in  robustness  theory  are 
applied  resulting  in  a  frequency  domain  condition  which 
guarantees  staoiiity.  The  stability  condition  is  applied  to  the 
design  of  a  shape/vibration  LQG  controller  using  position 


measurements . 


introduction 


I . 


Linea r- Quad ra t i c - Gaus s  i  an  (LQG)  control  theory  has  o een 
proposed  (for  example,  Refs.  1-4)  as  a  methoaolgv  for  designing 
control  systems  for  large  space  structurer  (LSSj.  (LQG  control 
is  taken  here  to  mean  linear  observer  oaseo  feedback).  The  basic 
difficulty  in  controlling  an  LSS  is  that  the  structure  is  lightly 
damped  (.2  to  .5%)  and  flexible.  Thus,  the  LSS  has  an  infinite 
number  of  highly  oscillatory  modes.  An  LQG  control  (or  any 
realistic  control)  is  finite  and  must  be  based  on  a  reasonably 
sized  finite  representation  of  the  LSS.  Consequently,  the  effect 
of  the  unmodeled  (residual)  modes  can  cause  instability. 

This  paper  examines  the  stability  of  LQG  control  of  the 
LSS.  Previous  work  along  these  lines  (for  example,  Ref.  S  )  has 
resulted  in  stability  bounds,  based  on  Lyapunov  Functions,  which 
can  be  conservative,  and  which  are  dependent  on  a  state-space 
representation  of  the  residual  modal  system.  Unfortunately, 
knowledge  about  the  residual  modes  is  uncertain  ana  might  be  more 
readily  available  as  a  transfer  function  matrix.  The  main  result 
presented  here  gives  frequency  domain  conditions,  which  if 
satisfied,  guaranty  asymptotic  stability. 

It  is  shown  that  the  effect  of  the  resid--al  modes  and  the 
observer  are  equivalent  to  placing  a  finite  gain,  linear 
perturbation  in  the  control  loop  of  an  LQ  (linear  quardratic) 
full  state  feedback  control.  Recent  results  [6-9]  in  robustness 
theory  for  multivariable  systems  show  that  there  are  frequency 
domain  conditions  which  guaranty  stability  of  control  systems  in 
the  presence  of  multiplicative  perturbations.  The  frequency 
domain  conditions  can  then  be  applied  specifically  for  LQG  modal 
control  of  the  LSS  and  can  be  used  for  design.  The  residual 
system,  i.e.  the  total  effect  of  spillover,  is  represented  as  a 
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transfer  function  matrix.  This  form  more  easily  allows  for  a 
limited  knowledge  of  the  residual  mooes.  For  example,  knowledge 
of  bounds  on  the  magnitude  of  the  residual  system  response  at 
different  frequencies  is  sufficient  to  test  for  staoility. 

The  paper  is  organized  so  that  Sections  II  and  III  state  tne 
basic  problem.  Section  IV  outlines  the  general  form  of  the 
frequency  domain  robustness  test  for  the  LQS  modal  control. 
Section  V  modifies  the  test  for  LQG  plus  feedtnrough  and  Section 
VI  shows  an  application,  in  general  form,  to  a  shape/vibrat  ion 
control  using  position  sensors. 
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II.  MODAL  REPRESENTATION 

Consider  a  finite  element  model  (FE.MJ  of  a  iigntly  lampec 
flexible  structure  given  oy , 


q  =  :- 


where 


are  -.-vectors  of  modal  coordinates  and 


generalised  displacement  coordinates  respectively;  u  is  an 
m-vector  of  point  control  forces;  is  an  x  d  1  agna  1  matrix  of 

modal  frequencies  ~  ^  ;  -  <  0  is  the  damping  ratio  {'■>  ~ 

.005  );  and  :  is  an  ■  x  i  matrix  whose  columns  are  the 

(approximate)  structural  mode  shapes.  Measurements  of 
displacment  and  velocity  are  produced  by  point  sensors  on  the 
structure.  Thus, 


+  Cvq 

V"  * 


(2.2) 


where  y  is  an  m-vector  of  measured  displacements  and 
velocities.  In  general  the  dimension  of  the  FEM  is  very  large 
and  equation  (2.1)  is  usually  decomposed  (by  an  appropriate 
ordering  of  modes)  into  P-primary  modes  and  R-residual  modes, 


i  .  e  .  , 


•’p  *  2  tp.ipPp  ♦  “pnp  3  Jp:U 


*■  ’R"RrR  ‘"R^R  *  :R  '  u 

y  =  CdVp  *  CvVp  +  CdVr  +  CvVr 


(2.3) 


w  n  e  r  e 


ana 


reoresen 


5 . Jua .  noaa. 


P 


;oorJna:e  vectors  respec 


tne  pri.iurv  ana 
v e  1  .  A  more  compact 


represenra: ; on 


xp  *  ApXp  *  ppij 

XR  =  ARxR  *  3Ra 

v  i-pxp  *  a x p 


3  ('p 

V 

X 

TO 

ii 

"R 

' R  ^  and 

0 

ip 

j 

A 

1 

r 

.  3 

P 

’  P 

.  T  _ 

1  - 

z 

l  p 

P‘ 

p  -i 

L 

P  J 

1 

- 

0 

lR 

0 

AR  = 

-> 

,  3r  = 

.T  - 

.  ‘-R 

_  1  - 
“■  ’  R” 

R  - 

1 

. 

'  R  '  . 

CP  =  ^  CD: P  C  v  *  P  ^  ’  CR  =  - CD*  R  C v*  R  ^ 


III.  LQG  MODAL  CONTROL 


The  LQG  modal  control  [  1  j  is  given  by, 

ULQG  =  'F  xp  (5.1) 

*  A  A 

*  ApXp  +  Bp  uLqG  +  K  (y-CpXp) 

where  (F,K)  are  selected  in  accordance  with  standard  LQG 
Theory  [  10  ]  based  upon  (2.4)  as  if  the  residual  modes  had  no 
effect.  Of  course  this  omission  is  only  an  artifact  to  select 
(F,K).  The  effect  of  the  resdual  modes  comes  through  the  terms 
Bpu  and  CpXp,  designated  [  1  ]  as  control  and  observation 
spillover ,  and  may  cause  degradation  of  closed-loop  performance, 
if  not  instability.  The  problem  to  be  addressed  in  this  paper  is 
to  find  conditions  which  guarantee  stability  of  (2.4)  using  the 
LQG  modal  control  (3.1). 
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be  r 1 ie 


e  r  3  <  p  -<p 

U  r  -  3  •  rX; 

l.  ;  r 

where  e  p  t  §  tie  pr  i.narv  noae  state  er'jr  estimate  a  no  ,• 
is  the  LQ  v  Linear  ,'uacratic;  full  state  p  r :  n  a  r  v  ncca  1  c  .nr  r  . 
which  c  ju  I  a  oe  used  if  the  rrmarv  moce  state  *ere 

available  tor  neasurement .  The  Lv,0  control  5..,  can  then  re 
expressed  as, 


Direct  . 

yield , 

ULQG  =  ULQ  '  Fe? 

^  p  3  ( A  p  *  KC  p ;  e  p  * 

:alculat ions  using  Laplace  transform  as 

>.3.3,  ana  1  - 

where 

uLQG(s)  =  Lls^  uLQis) 

Lis)  =  (I  *  c i s ) ) 

E  (  s  )  =  F  Jp (  5  )  K  ( I  *  C p  ip ,  s  J  k  ;  1  Hp  ■,  s 

3.3 

HR( s )  =  Cpi R( s ) 3R 

with 

ip(s)  3  (sip  -  Ap) 

}R(s)  =  ( s I R  -  Ar)-1 

i  3  o 

Equation  (3.5)  can  be  interpreted  to  mean  that 
control  (3.1)  of  the  modal  decomposed  system  ( 
to  a  multiplicative  perturbation  of  the  ideal) 

the  L  Q  G  moaai 

2.4)  is  eauivalen 

ted  LQ  modal 

control.  Figure  1  illustrates  the  block,  diagram  equivalence. 
Consequently,  stability  analysis  of  LQG  modal  control  can  be 
viewed  as  robustness  of  LQ  modal  control  to  a  perturbation 
represented  by  the  transfer  f  unc  t  ion  Lis)  in  the  input  c.hanne 

IV.  ROBUSTNESS  TEST  FuR  STABILITY 


Recent  studies  [6-9]  of  the  stability  and  rooustness  of 
multivariable  systems  has  shown  that  for  the  perturbed  system  of 
Figure  2  ,  if 

£  l  ♦  >7  i:G(>)  ,  ~£0  14.1; 

then  the  perturbed  system  is  stable.  The  quantities  £(  •  )  and 
7(0  are  the  minimum  ana  maximum  singular  values  respectively, 
of  the  matrix  argument.  equating  the  quantities  in  Figure  2  to 
Figure  1  gives, 

G(s)  =  FJpCsJBp  U.2J 

dG(  s )  =  Us)  -  I  »  -Cl  *  E'l(s)  )'1 

since  7  ( M "  ^  )  =  l/c(M)  for  any  matrix  M  then  [4-1)  with 
(4-2)  becomes, 

£  ( I  *  G’l(jii))  £(I  *  E‘l(j-)j  >1,  ->0  (4.3) 

Thus ,  it  inequality  (4.3)  is  satisfied  then  the  LQG  modal  control 
stabilizes  the  LSS. 

An  important  aspect  of  the  robustness  test  (4.3)  is  that  it 
is  a  frequency  domain  test  and  the  spillover  effect  is  contained 
in  E(s)  in  the  term, 

Hq(s)  *  CqiR(s)BR  (4.4) 

This  term  is  the  residual  input-output  transfer  function  matrix. 
In  a  practical  sense  knowledge  of  bounds  on  Hp(s)  may  be  more 
readily  available  then  the  state  representation  (aR,  3R> 

Cp)  particularly  for  the  high  modal  frequencies. 

A  more  conservative  test  than  (4.3)  can  also  be  developed. 

It  is  shown  in  [11]  that  for  LQ  control, 
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consequently  ,  -l  .  5  j  is  surely  satisnea  i  f 
"vEi.  )  j  _  1  ,  ..  >  0 

See  Appendix  A  for  an  alternate  derivation  oasea  on  Ret. 


A  .  5  / 

.  A  .  o  , 

3. 


V.  LQG  WITH  DIRECT  CUTPUT  FEEDBACK 

The  preceding  analysis  is  also  compatible  for  LQG  plus 
direct  output  feedback.  That  is,  consider  the  control, 

u=^-LXp-Dy  i  s  •  1  j 

using  (2. A],  (3.5),  and  (5.1)  results  in, 

u  =  -  L(.Xp  ♦  ep)  -  D  ( CpXp  CRxRj 

=  -  i  L  +  DC )  Xp  -  Lep  -  DCRxR 

=  ULQ  "  ^ep  '  DCR  xr  (5.2j 

w  hers 

ULQ  3  *  Fxp>  F  3  L  +  DC  15.3) 

Thus,  the  resiaual  and  error  dynamics  are, 

*R  =  aRD  xR  ■  BRLeP  *  3ruLq  (5.A) 

ep  =  (Ap  -KCpJep  ♦  KCR  XR 

where 

ARD  3  AR  ‘  BRDCR  (5.5) 

According  to  Canavin  [2]  the  feedthrough  gain  D  is  used 
primarily  to  influence  the  residual  characteristics  AR^ 

The  perturbation  transfer  function  now  becomes  (See  Appendix  B), 
LD(s)  •  (I  -DHRD(s))  (I  -  Q(s)Hrd(s)  T1  1 5  . 6  ) 

where 
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UR  'rd1 s  j  3r 

1  S I R  ‘  ARD ^  1  3  ‘  ' 

L  1  c  5  )  K  l  I  *  C  :  l  s  )  K  J  " 1 
?  p  p 

If  E  t-j  C  s  )  is  defined  so  chat, 

Lq(s)=(I*Eq(sJJ*  (5.3; 

then  using  (25)  and  v.  2  ~ )  results  in, 

Ed(s)  =  ( Q(  s )  -  Q)  Grd(sKI  -  DGrd1s))'1  (5.9; 

Following  (4.5),  if 

£,(  I  G  ^(.j-j))  £_  C I  ♦  1,  *._>  0  (5.10; 

then  LQG  plus  direct  feedthrough  stabilizes  the  system  (2.4). 
Similarly,  following  (4.6),  if 

■£ED(j-)J  1  1  ,  -io  (5.11) 

the  LSS  is  also  stable. 

VI.  APPLICATION:  SHAPE/VIBRATION  CONTROL 

The  robustness  test  (4.6)  will  now  be  applied  to  an  LQG 
shape  and  vibration  controller.  Assume  that  there  are  p 
colocated  force  actuators  and  position  sensors  on  the  structure. 
It  is  desired  to  control  p  primary  modes  using  an  LQG  approach 
which  neglects  (for  the  design)  the  residual  modes.  Thus,  the 
primary  mode  system  matrices  (Ap,  Bp,  Cp )  are 


hrd!s) 


rRD(  S  ’ 


Q(s)  = 
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After  nuch  algebraic  manipulation  :  anc  using  me  race  r.na 
diagonal  matrices  commute;  the  terms  in  ,3.5,  oectme, 


;o'^;K  =  V  3  >;'5;  ,H5':o 


O  t) 


S  ;  .\  = 


D  *  I  S  ;  ,  S  I  7  , 


The  matrices  Dts,),  Vis),  -vis)  are  ail  pxp  diagonal  a 
given  by, 


lKs)  =  5*’ 


P  ■'? 


Wis)  =  2s(;0.10(.T- 


■  Z) 


1 i 


r  -  _  - .  «  ■)  , -  a-  -  "  ■  '  “ 

c  "  0  ""  o  "?  "*  ’  C’  0"  C"  0"  p 


v(  s )  =  s'  l„  -  s  2*:. 


'  p  0  "0  '  0 

Using  16.6)  and  (6.7)  with  (3.5)  gives 


*  1  n  *  1  c  I'n  =■  >  -  *  1  1  -  1  n  *  ^  i  c  i ;  c  ' 


Eis)  =  :"d  i-(s)vus):„  i:‘D  hsjV.s  J 


o  o 


?  '  R 


m0  1  s 


T'V 1 ( s) W( S ) V 1 ( S )D( s ):;T  Ha ( 5 ) 

?  r  A 


■W(s)V*1(s)r*THRis) 


C  rive  s )w  ~ 1 C  s)r  ^ ) ' 1  hr(s) 


? . . . .  p ' 

Consequently,  (  a  .  6  )  is  surely  satisfied  if 


r  ( H9  (  jij ) )  <  r  2  ( 7  „  )cr  (V(  jo,  )W  Li  ju  ) )  ,  w  >  0 


w-1- 


*R"  '  '*  p- 

From  (6.3)  ana  (6.9)  the  diagonal  elements  of  V( s )  and 


are 
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whe  re 


V  ■  Is)  *  s'  ♦  s  2  -  ~ 

J  3  '’0“0  o 


Vs)  ■  a t s  -  lij 


;  I 


i  =  1,2 . p  '0.15; 


ai  =  2(  Vo.  (jc.  *  “p.}  +  ^cJc-Lo.  '“P  3 

7  7  ?  ?  ? 

^  i  ^  M  C  "n  -1  p  ^  ^  -  c  0"  C  “0  ~  p 

i  *i  i  r  l  i  l  i 


v.  6  .  1 4  2 


Thus , 


r  (Hr(  jcj)  )  <_ 


min 
i-1. . . 


(  6  . 1 5  J 


for  all  ~  >  0,  guarantees  stability. 

The  left-hand  side  of  (6.15)  must  be  bounded  for  all 
frequencies.  This  requires  that  the  LSS  have  some  natural 
structural  damping.  Assuming  damping  to  exist,  there  is 
guaranteed  stability  of  the  LQG  control  if  the  right-hand  sice  is 
unbounded  for  all  frequencies.  This  condition  would  require  all 
a^  and  b^  to  vanish.  One  way  to  achieve  this  is  to 
eliminate  the  observer  (K  =  0)  and  use  rate  feedback  only 
(F^=o,  F2  =r ‘1Q,  Q=QT  >0)  (see  Refs.  12  and  13).  This  is 
a  limiting  case  because  there  is  no  observer  and  no  shape 
control,  however,  (6.15)  is  validated  with  respect  to  a 
well-known  result. 

In  general,  it  does  appear  that  (6.15)  can  be  almost  always 
satisfied  by  some  selection  of  controller/observer  gains.  The 
gains  are  thus  constrained  and  there  is  a  trade-off  between  the 
desired  behavior  of  the  primary  modes  and  robustness  to  spillover 
( res idual  ef fects) . 
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Another  point  worth  noting  is  that  the  right-hand  sue  of 
(6.15)  is  quite  accurately  known  since  it  is  created  oy  aesign. 
The  left-hand  side,  which  represents  the  spillover,  nay  not  be 
entirely  known  for  all  frequencies.  However,  knowledge  of  bounds 
of  the  residual  transfer  function  at  different  frequency  oands 
satisfies  the  information  required  in  (6.15).  For  example,  given 
the  problem  being  examined,  the  residual  mode  system  matrices 


(A 


R>  BR 


cr) 


are 


AR 


0 

-> 

*R 

3  =  r^-  = 

0 

.  *nR 

4’R"R 

> 

°R  UR 

.  :R  . 

Consequent ly , 


Hr(  s) 


cr1>r(s)  br 


:J(s2Ir  ♦  s  2.; 


R“  r 


V1  - 

■R;  •  i 


Thus 


'(Hr(  j  J)  < 


r2( :R) 


min 

i  =  l .  .  .R 


A  2 


“R. 


2)2  +  4.2  2  2 

'  1  -R*  Ri- 


(6.16) 


Combining  (6.16)  with  (6.15)  gives. 


o(:R) 


min 

i  =  l .  .  .R 


min 

i  =  1 .  .  .  p 


t  2  2)'  .,2  2  2 

UR.  -0  '  -R^R.^ 


f  2  2  ,  ,  2  2  2 


,2  2  ,2 
ai  +  bi 


(6.17) 
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Use  of  an  expression  as  specific  as  . o  .  ! '  ■  or  .:  )  u  r  s  e  requires 
chat  (6.16)  is  accurate,  particular!*'  ->  t  cue  nign  frequencies. 

VII.  CONCLUSION 

The  stability  of  LQG  modal  control  tor  an  LSS  nas  oeen 
examined.  Results  from  recent  researcn  in  rooustness  tneor-', 
gives  rise  to  frequency  domain  tests  tor  stability.  These  test 
can  be  used  even  in  the  face  of  uncertainty  aoout  t ne  re s : uua . 
modes,  particularly  the  higher  order  mooes. 

Future  research  using  the  methodology  presentee  here  sr-.oul 
examine  (1)  the  effect  of  feedthrough  (e.g.  colocateo  active 
dampers)  on  the  LQG  control,  (2)  the  effect  of  either  sensor  3 r 
controller  prefilters,  and  (3)  the  effect  of  sensor  ana 
controller  dynamics.  It  also  appears  that  the  question  of 
decentralized  control  of  LSS  could  be  examined  m  tms  wav. 
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APPENDIX  A 

ALTERNATE  DERIVATION  JF  ROBUSTNESS  TEST 


Following  [9]  an  LQ  control  will  be  stability  rooust  for 
finite-gain,  linear  input  multiplicative  perturoat ions ,  with 
transfer  function  L(s)  if 

U  iu  )  *  LT(.  -  j*  )  _>  I  ,  ^  0  ( A .  1 ) 

From  (3.5) 

UsJ^U+Eis))1  v  A .  2  ) 

Using  (A. 2),  ( A . 1 )  becomes, 

( I  *ET(  -  i  ~)  )  [  LC  jw  )  *LTC  -  j*  j  1T(I+E  t-jw)  >  (l*ETl-  jjJ  KdElj^J  ) 

Thus 

21  -  ET(-jw)  =  EC  j-i)  >  I -ETL  -  J  «■  ECjy)  ♦  E T  (  -  j  w  J  E  C  j-- J 

and  finally 

ET(-j^)  E(j^)  _<  I  (A.  3) 

which  is  implied  by 

TCEC  ju.)  )  <1,  ,  >  0  (A.  4 ) 


163 


APPENDIX  3 

LQG  PLUS  DIRECT  FEEDTHROUGH 

To  develop  (5.6)  requires  the  transfer  function  between 
xr  and  u^q  and  ep  and  xR.  Direct  calculation  from 
(5.4)  gives , 

XR  “  * 1  *  ®rQ' 3 )Cr  J  1  uLQ  t3' 

ep  -  (Jp1(s)  *  KV’1  kcr  xr  u- 

where  Q(s)  and  JrqCs)  are  given  by  (5.7).  From  (5.2) 

ULQG  =  ULQ  ‘  L  ep  '  DCR  xR 
3  ULQ  '  (Q(s)>D)C  xR 

-  [I-(Q(s)-D)CR(I*?RD(s)B.lQ(s)CR)'1tRD(sj3R]uLi, 
=  Lq(s)  uLq  (B. 

Using  the  matrix  inversion  lemma  on  LqIs)  gives, 

Ld(s)  3  I  -  ( Q  +  D  )  CR  [  I  -  5  RQBR  ( I  +  Q  Cr^Br)  Q  Cr^rd3R 

■  I-CQ»DHHrd-HrdCI»Q  Hrd)'1Q  Hrd) 

■  I-CQ*D)HRD(I.Q  hrd> " 1 

■  (I-D  HRD)(I.Q  Hr0)-‘  IB. 
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